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Abstract 



The masses and decays of the scalar D*o(2317) and axial- vector D*]^(2460) charmed 
strange mesons are calculated consistently in the hadrogenesis conjecture. These 
mesons decay either strongly into the isospin-violating vr'^Ds and vr^D* channels or 
electromagnetically. They are generated by coupled-channel dynamics based on the 
leading order chiral Lagrangian. The effect of chiral corrections to chiral order is 
investigated. We show that taking into account large-Nc relations to determine the 
strength of these correction terms implies a measurable signal for an exotic axial- 
vector state in the r] D* invariant mass distribution. The one-loop contribution to 
the electromagnetic decay amplitudes of scalar and axial- vector states is calculated. 
The Lagrangian describing electromagnetic interactions is obtained by gauging the 
chiral Lagrangian for hadronic interactions and adding gauge-invariant correction 
terms to chiral order Q^. In addition the role of light vector meson degrees of free- 
dom is explored. We confront our results with measured branching ratios. Once 
the light vector mesons are included, a natural explanation of all radiative decay 
parameters is achieved. 
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1 Introduction 



The observation [Tp] of two narrow, positive parity, charmed strange mesons at 
masses lower than expected in quark models [31H] may provide new insight into 
the way hadrons are generated. The properties of these two mesons, the scalar 
D*q(2317)^ and the axial-vector D*]^(2460)^, appear indeed sensitive to the 
degrees of freedom building up hadronic excitations and to strong interaction 
symmetries [511^ . It is therefore of interest to compute the properties of these 
states in effective field theories involving hadronic degrees of freedom and 
constrained by specific symmetries. 

The importance of symmetries in constructing effective actions describing 
heavy-light mesons was emphasized a long time ago. As early as 1993, Nowak, 
Rho and Zahed [7j derived an effective action combining chiral and heavy- 
quark symmetries and predicted that the pseudoscalar (0~) and vector (1~) 
ground states should have chiral partners as a consequence of the spontaneous 
breaking of chiral symmetry. The splitting between the two sets of states was 
found to be rather small (of the order of the constituent quark mass). It was 
also noted that the large Nc limit appeared compatible with the heavy-quark 
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limit in the heavy-light meson sector. Later work by the same authors [8] 
addressed specifically the observed D*q(2317) and D*^(2460) mesons. Similar 
considerations linking chiral and heavy-quark symmetries to parity-doubled 
heavy-light mesons systems were published at the same time by Bardeen and 
Hill |9] . These papers rely on the chiral quark model which predicts the heavy- 
light 0"^ and l"*" states as chiral partners of the heavy-light 0~ and 1~ antitriplet 
states. There is no prediction of that kind for non-linear realizations of the 
chiral SU(3) group. 

The discovery of the D*o(2317) and D*i(2460) mesons in 2003 and the diffi- 
culty of conventional approaches to reproduce their masses motivated a large 
number of theoretical studies involving different descriptions based on multi- 
quark states, molecular pictures or dynamical generation of resonances (see 
[511^ for a broad set of references) . 

In order to get more insight into the nature of the D*q(2317) and D*]^(2460) 
mesons, it was pointed out by Mehen and Springer [10] that the electromag- 
netic decays of these states could be most helpful in distinguishing among 
models. This property was illustrated by a specific comparison of leading or- 
der heavy-hadron chiral perturbation theory predictions for these decays to a 
corresponding calculation in the molecular picture. 

The purpose of our paper is to study consistently the masses, electromagnetic 
and strong decays of the D*q(2317) and D*^(2460) mesons in the hadrogen- 
esis conjecture where earlier work [TT]|12fl3] showed that such states can be 
produced at their observed masses. This approach for heavy-light mesons ex- 
ploits both heavy-quark and spontaneously broken chiral symmetries and ge- 
nerates D^-mesons through relativistic coupled-channel dynamics. Goldstone 
bosons are scattered off heavy-light pseudoscalar and vector D-meson ground 
states. For the D*q(2317)"'" meson, the calculation involves the rjD'^, K^D^ 
and K^D^ channels coupled further to the tt^D'^ channel through an isospin- 
mixing parameter. Analogously we consider for the D*^(2460)"'" meson the 
r]D*~^, K^D*^, K^D*^ and 7r°_D*+ channels. The strong, parity-violating de- 
cays of the D*q(2317) and D*^(2460) mesons give the main part of their total 
widths and reflect both the coupled-channel dynamics and the treatment of 
isospin-mixing effects. The electromagnetic widths are very sensitive to the 
details of the coupled- channel dynamics and to the vector degrees of freedom 
treated explicitly in the interaction Lagrangian. 

A series of papers have addressed specific issues related to the contents of 
our work. The D*q(2317) and 0*^^(2460) mesons were dynamically generated 
based on the leading order heavy chiral Lagrangian by Guo et al. [T^IITS] . They 
also computed the strong decays of the D*q(2317) and D*;^(2460) mesons. As 
this calculation is quite comparable to our first step leading order deriva- 
tion, we will compare our results for the strong widths and comment on the 
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differences. Open-charm meson systems and their decays were also studied 
from a Lagrangian based on the SU(4) flavour symmetry by Gamermann et 
al. [TB]|17fl8] . The starting SU(4) degeneracy in this approach is quite dif- 
ferent from the light-quark chiral symmetry and charmed quark heavy-quark 
symmetry underlying our work. The SU(4) symmetry is not an approximate 
symmetry in the presence of a heavy quark and needs to be largely broken 
by phenomenological interactions [16]. It is therefore expected that the prop- 
erties of heavy-light mesons will differ significantly in both approaches even 
though specific effects in radiative decays can be similar. We will illustrate this 
point in the discussion of the electromagnetic decay of the D*q(2317) state. 
The D*q(2317) and D*^(2460) mesons were interpreted as hadronic molecules 
by Faessler et al. p^]|20] . It is assumed that the D*q(2317) is a strong bound 
state of K and D mesons while the D*]^(2460) is treated as a bound state of 
K* and D mesons. The strong and radiative decays were calculated using a 
phenomenological Lagrangian. While this phenomenological approach is quite 
far from our description, we find some common features in the strong decay 
properties. We will point them out. 

We summarize now the experimental data available on the charmed strange 
mesons. We display in Fig. 1 the D^-meson spectrum as presently known 
[2T|22] . We also indicate the DK threshold whose closeness to the D*q(2317) 
and D*^(2460) mesons influences the dynamics of these states. The spin and 
parity of the D^-mesons are well-established for the ground state and for the 
D*]^(2460)^. The spin and parity of the other states need confirmation. We 
have quoted their most probable values. 



The D*(2112)=^ and D*o(2317)=^ mesons lie below the DK threshold and are 
therefore expected to be very narrow states. They can decay either electro- 
magnetically or into the isospin-violating vr^ channel. 

The D*(2112) has a width F < 1.9 MeV and decays dominantly by a radiative 
transition to the scalar ground state with a probability of (94.2 ± 0.7) % [2T]. 
Its decay probability to the Dgir^ channel is therefore (5.8 ± 0.7) %. 

The most stringent upper limit obtained for the D*q(2317) width is F < 
3.8 MeV [23]. The D*q(2317) was first observed through its decay into the 
0^(1968) 7r° channel [T]. Its radiative decay to the 0^(1968) has never been 
seen. The upper limits available on the ratio of the radiative to pionic decay 
widths of the D*o(2317) to the ground state and to the D*(2112) are [21] 



si 
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(94.2 ± 0.7) % 



2573 MeV 1=0 J"=2^ 
2536 MeV 1=0 j'^ = 1 + 



2460 MeV 1=0 j"= 1 + 



DK threshold 

2317 MeV 1=0 J"=0^ 



21 12 MeV 1=0 j"=l 



1968 MeV 1=0 J'^ = 



Fig. 1. The D^-meson spectrum with the most probable spin-parity assignments 
together with the radiative transitions on which there is experimental information 
|21|22j . The arrows without number indicate that there are experimental constraints 
on these decays but no absolute branching ratios. The dashed line shows the DK 
threshold. 



r [L):o(2317) ^ L>,(1968) 7rO] 



< 0.059. 



(2) 



The D*^(2460) meson is located above the DK threshold but appears never- 
theless very narrow: its total width was found to be less than 3.5 MeV |23j . 
Constraints on its radiative decays to the Ds(1968), to the D*(2112) and to 
the D*o(2317) are as follows [2T], 



r[D:,(2460) -> D,(1968)7] 
r[D*i(2460) ^ D*(2112)7r0] 



0.31 ±0.06, 



(3) 



r[D:i(2460) ^ (2112) 7] 
r[L)*i(2460) ^ L)*(2112)7r0] 



< 0.16, 



(4) 



r[J:,(2460) ^ /^:o(2317)7] 
r[L)*i(2460) ^ L)*(2112)7r0] 



< 0.22. 



(5) 
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An absolute measurement of the decay probability of the D*]^(2460) to the 
Ds{1968)~ 7 channel gave recently (16 ± 7)% [22]. The same data also provided 
the branching fraction B{D*^{2AQ0y D;(2112)~7r°) = (56 ± 22)%. The 
ratio is in agreement with the value quoted in 

The D*i(2536) and D*2(2573) widths were recently measured to be (1.03±0.017) 
MeV and (27.1 ±6.2) MeV respectively [2^||25] . The radiative decay widths of 
these states are not known. 



We restrict our calculations to the radiative and isospin-violating vr^ decays 
of the D*q(2317)^ and D*]^(2460)"'" states on which there are fragmentary but 
significant data. We note that the presence of a 2"*" state about 100 MeV 
above the D*^(2460)^ state is an indication of physics outside our scheme 
around that energy scale which should guide future investigations of higher- 
lying states. In view of the uncertainties in the measured widths of the D*q and 
D*i states, our main concern will be to check the ability of the hadrogenesis 
conjecture to provide a consistent picture of the main features of the decay 
scheme. We aim at identifying the important contributions to the dynamics 
of these decays, determining characteristic ranges for the parameters involved 
and making predictions able to test further the structure of the D^-mesons 
and the specific conjecture on which this work relies. 

This paper is organized as follows. Section 2 is devoted to the generation of 
the D*q(2317)+ and D*^(2460)+ mesons and to the calculation of their strong 
isospin-violating decay widths. We present first a calculation based on the lead- 
ing order chiral Lagrangian in which massive vector particles are described in 
terms of antisymmetric tensor fields. We introduce subsequently chiral correc- 
tion terms to chiral order to take into account s- and u-channel D-meson 
exchange processes and local two-body counter terms. Section 3 deals with the 
coupling of the electromagnetic field to the hadrons. In a first step we gauge 
the hadronic interactions introduced in Section 2. We add gauge-invariant in- 
teraction terms of chiral order Q^. We consider also interaction vertices probed 
when including the light vector mesons as explicit degrees of freedom. We com- 
ment on the values of the parameters associated with these terms in relation 
to QCD symmetries and discuss the renormalization of the ultraviolet singu- 
larities. The explicit expressions of the electromagnetic decays are derived in 
Section 4 for the scalar state D*q(2317) and in Section 5 for the axial-vector 
state D*^(2460). Our numerical results are presented in Section 6 and com- 
pared to the available data. We discuss the role of the different contributions 
and the constraints expected on the range of values for the coupling constants 
of specific interaction terms. We conclude in Section 7. We relegate lengthy 
derivations in seven appendices (A-G). 
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2 Generation and strong isospin-violating decays of molecules 

This section deals with the generation of the scalar -D*o(2317) and isovec- 
tor D*i{24Q0) mesons in the coupled-channel framework of p!2l[T3] . This des- 
cription is based on the scattering of Goldstone bosons off heavy-light 0~ and 
1~ mesons respectively. Isospin-breaking effects arise from the difference be- 
tween the up and down quark masses which leads to isospin-violating strong 
decay amplitudes, D*q{2317) Ds and D*i(2460) — > D*. These isospin- 
breaking effects were not included in earlier work based on the assumption of 
perfect isospin symmetry [T2IT3] . We consider scalar and axial- vector states 
successively. In the latter case, we reformulate the derivation of [T2IT3] in terms 
of massive 1~ fields represented by antisymmetric tensors. This particular de- 
velopment is needed later to arrive at expressions which are gauge-invariant in 
a transparent manner. For both states, we consider first the hadronic interac- 
tions resulting from the leading order chiral Lagrangian and treat afterwards 
chiral corrections to chiral order Q^. 

We emphasize the importance of studying scalar and axial-vector mesons in 
the open-charm sector on equal footing. The properties of spin and spin 1 
heavy-light mesons are indeed closely related by the heavy-quark symmetry of 
QCD P^27f28f29ll3U] . Even though the charm quark mass is much larger that 
the light (u,d,s) quark masses, the limit in which the mass of the charm quark 
goes to infinity is an approximation which may require significant corrections 
[3T|32] . Rather than applying a formalism where scalar and vector fields are 
fully degenerate and grouped together in one field as implied by exact heavy- 
quark symmetry, we use separate scalar and vector D-meson fields to allow for 
the observed mass difference between spin multiplets. 



2.1 The scalar state D*q{2317) 

The open-charm D*q(2317) state has been shown to be dynamically generated 
as a direct consequence of the leading order chiral Lagrangian density [T2IT3] . 




where $ and D are the pseudoscalar octet and triplet fields. We use the no- 
tation D = D"^ . In the particle representation the Goldstone and ground state 
open-charm meson fields are 
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TT 



V 



V2 7r- 
y/2K 



(7) 



The Weinberg- Tomozawa term in ([6]), which is proportional to /~^, is obtained 
by chirally gauging the kinetic term of the D-mesons. It is of chiral order Q^^ 
as it involves a single derivative of the light Goldstone fields. The parameter 
/ ~ /tt = 92.4 MeV in ([6]) is the octet meson decay constant. It defines 
the scale of chiral symmetry breaking and is approximatively known from 
the weak decay of the charged pions. A precise determination of / requires a 
chiral SU(3) extrapolation of some data set. In [33] the value / ~ 90 MeV was 
obtained from a detailed study of pion- and kaon-nucleon scattering data. We 
will use / = 90 MeV throughout this work unless specified otherwise. This 
parameter determines the leading s-wave interaction of the Goldstone bosons 
with the open-charm meson fields. We note that our approach is consistent 
with heavy-hadron chiral perturbation theory [10] but at variance with the 
work of [T6|17II18] where the large breaking of the assumed SU(4) flavour 
symmetry leads to different interactions. The ground state scalar D-meson 
mass matrix is denoted by Mq-. The mass term of the Goldstone bosons is 
proportional to the quark-mass matrix 



^ 



Xo = 25o 





m. 



2 m 



K 



ml)Xs. (8) 



At leading order the latter can be expressed in terms of the pion and kaon 
masses as indicated in ([8]). 

If we admit isospin-breaking effects, i.e. m„ ^ rriii, there is a term in ([6]) 
proportional to {mu—md) 7^^ rj, inducing vr'' r] mixing. A unitary transformation 
is required such that the transformed fields vf" and rj defined by 



7r° = 7r° cos e — fj sin e , 



7] = 71^ sin e + ?7 cos e 



(9) 



decouple. The Lagrangian density ([6]), when written in terms of the new fields, 
does not show a. fj term if and only if 



sin(2e) 



rrid - rriu 



cos(2e) ' 2ms—mu—md 
According to [M] the ratio of quark masses relevant in (ITOl) takes the value 



(10) 
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which imphes the mixing angle 



e = 0.010 ±0.001. (12) 

Heavy-hght meson resonances with quantum numbers = 0+ manifest them- 
selves as poles in the s-wave scattering amplitude T{s). We consider the four 
isospin states {KD,I\, {7iDs,l\ and {riDs,0\. In the presence of isospin mi- 
xing all channels couple. The mixing of the two isospin sectors is of order e. 
Using the phase convention of [T2IT3] . we define the four states 



(ll 




= cose (tt^ Dg, 1 


+ sine {r] Dg, 


(2| 


= {V Dt\ 


= cose {r] Dg, - 


- sin e (tt^ Dg, 1 


(3| 


= (/i°Z}+| 


= +^((irAO| 


~{KD,1\), 


(4| 


= {K+D'\ 


= -j=^{{KD,0\ 


+ {KD,1\). 



(13) 



The Weinberg- Tomozawa interaction (jSj) implies a scattering amplitude of the 
simple form p!2l[T3] 



T{s)=[l-V{s)J{s)] W{s). (14) 
The matrix of loop functions, J{s), is diagonal and given by [T2p!3] 



As) - 




m 



(15, 



where ^/s = y -|- + y -|- p^^^ and Pcm is the center of mass momen- 
tum of the scattering particles. Each diagonal element depends on the masses 
of the light (Goldstone) and heavy (open-charm) mesons m and M defining 
the coupled-channel state. The matching scale /im in f[T5]) should be identi- 
fied with the ground state mass of the D^-meson, i.e. fiM — 1968 MeV |33j . 
For such a value, s- and u-channel unitarized scattering amplitudes may be 
smoothly matched around the matching scale so as to define a full scattering 
amplitude that is crossing symmetric by construction. The determination of 
the matching scale and the induced approximate crossing symmetry parallels 
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the derivation published for kaon-nucleon scattering in [33]. The minimal crit- 
ical point needed to open the matching window for vrD., or KD scattering is 



typically of the order of + M|,^ ~ \l'raK^ + M|, ^ Mj:,^. One may vary 
the matching scale slightly around its natural value. As shown in [11] the re- 
sulting effects are small for reasonable variations. A large variation is excluded 
since it would make the matching of u- and s-channel unitarized amplitudes 
possible only at the price of introducing a strong discontinuity, which is at 
odds with causality. 

The coupled-channel interaction kernel Vij(s) in f|T^ is determined by the 
leading order chiral SU(3) Lagrangian ([6]) to be [12] . 



T/(0+) , . _ ^ 

VwT,ij{S) ~ gp 



X 



where (m, M) and (m, M) are the masses of initial and final mesons and 
the indices i and j refer to the states defined in f|T3l) . The 4x4 matrix Cjj, 
whose elements characterize the interaction strength in a given channel, can be 
expressed in terms of the mixing angle e and the isospin zero C\f and isospm 

one Clf' coupling matrices of [12]. For the channels of positive strangeness 
considered in ffT3l). we have 



Cn = Cfi sin^ e + C^^ cos^ e C12 = (Cif - C[l^) sin e cos e 

C22 = Cif cos^ e + C[^^ sin^ e ^13 = ^ (cg^ sin e - cos e ^(2^ 

C23 = ^ ( cos e eg' + sin e C^) C33 = 1 (cff + C^) 

Ci4 = ^ (Cf? sin e + eg) cos e) C24 = ^ (cg^ cos e - Cg' sin e 

C[f = 2 = 2 eg) C[f = V3 Cfi = eg) = eg) = . (17) 



Given the coupled-channel scattering amplitude ( fMl) with the effective inter- 
action ( fT6l) . it is straightforward to determine the mass and width of possible 
resonances. The Weinberg- Tomozawa interaction is strongly attractive in the 
isospin strangeness (I,S) = (0,1) sector, leading to the formation of a scalar 
resonance of mass Mq+. The latter manifests itself as a pole in the scattering 
amplitude which factorizes close to the pole at s = Mg+, 



J- i 
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with the couphng constants gi and the width ro+. We do not consider the 
possibihty of background terms as we are deahng with a narrow resonance. 

The numerical resuhs obtained at leading order should be viewed as quali- 
tative. At the end of this section more quantitative results including chiral 
corrections will be presented. 

To reproduce the empirical mass of 2317.6 MeV at leading order for the O"*" 
state requires using an effective parameter /e// = 95.5 MeV in f|T6|) . Taking 
e = 0.010, the coupling constants are 



QnDs - 1-95, g^oo, ^ 0.056, 

gKOD+ - 2.25 , gK+D^' - -2.25 . (19) 

Isospin-breaking effects in the KD coupling constants are found to be negligi- 
ble, i.e. gKOD+ — ~9k+d° holds quite accurately. It is therefore meaningful to 
work with the isospin coupling constants 



= V2gK^o+ ^ 3.18 , 4°+] = g.,D. ~ 1.95 . (20) 

The flavour SU(3) limit suggests ^/^g^^o} = g^Ku , a result quite compatible 
with the values given in f|T9| [201) . 

The strong width of the D*q(2317) is related to g.„oj^^ by 



TDi,[2m^.^D. = \g^onf'^^ |(7^o,5j2 23.7MeV ~ 76keV, (21) 

where pcm is the pion momentum in the center of mass frame. Our value of 76 
keV is almost an order of magnitude larger than the value given in [13] based 
on the same interaction. 

To trace the origin of this difference, it is useful to understand the physics 
underlying gnOOs- At linear order in isospin breaking, it is the sum of two 
terms. 



9.^D. = ~e9^KD^+egXK (22) 

The contribution proportional to g^j^^ is unambiguously determined by the 
angle e, as a direct consequence of the n^r] mixing defined in (Q, and linked 
to the finite value of (m^ — m^). If it were the only contribution to g^oDsJ 
the corresponding strong width Tf)*^(^23n)^TTODs would be ~ 9 keV, a value 
compatible with the 8.7 keV obtained in [13] based on n^r] mixing. 
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The contribution proportional to g\^j^ appears at the same level and is not 
included in |T3] . It is induced by the mass difference between the neutral and 
charged kaons and between the neutral and charged D-mesons. Like the ti^t] 
mixing phenomenon, it couples the two isospin sectors and is related to the 
finite value of (m^j — m^). It is clearly not there if one assumes mx+ = mj^o 
and M£,+ = M^o. This effect in chiral coupled-channel dynamics is analo- 
gous to the mixing phenomenon recently pointed out in the molecular picture 
p!9|20] where it is induced by the exchange of vector K*- and D*-mesons. The 
parameter e in ( l22l) must have the form 

6 = 61 + 62, 61 ~ (m^+ - mj^o) , 62 ^ {Mu+ - Mdo) , (23) 

where the proportionality factors depend on the details of the coupled-channel 
dynamics. At leading order we predict e ~ 0.012 with ei ~ 0.004 and 62 — 
0.008. The term eg^j^~l^ in g-^oij^ is roughly twice larger than eg^j^^ and appears 
dominated by the isospin-breaking induced by M£,+ 7^ M^o. These results 
agree qualitatively with the findings of [T91I20] . 




We mention two sources of uncertainty in the value derived above for the 
strong width of the D*q(2317). First we use the physical masses of the neutral 
and charged kaons and D-meson including the electromagnetic contribution 
to these masses. The latter should in principle be generated by the coupling of 
the hadronic Lagrangian to the electromagnetic field. The splitting of hadronic 
and electromagnetic interactions is a nontrivial issue [35)136] which we do not 
address here. The use of physical masses in the unitarization loop function 
(1151) is in line with the scheme developed in [33]. Loop corrections implied 
by photon-exchange processes enter the effective coupled-channel interaction 
V{s) in f|T^ at order Q^. Another source of uncertainty lies in the coupling 
constants gr^o^ and gxD- It was shown in [TH] that chiral correction terms 
of order modify these couplings and lead to the values grjOs — 3.7 and 
gxD — 3.7. The width of 76 keV quoted above is therefore most likely a lower 
limit as suggested by inserting the coupling constants of [13] into (!2T|) and 
(1221) . We will return to this issue in the final part of this section. 

2.2 The axial-vector state D*i(2460) 

We turn to the generation and strong decay of the axial- vector meson D*i{24Q0) 
This calculation involves the same procedure as followed for the scalar state 
except that we now build the axial- vector state by scattering Goldstone bosons 
off vector D-mesons. 

Since we aim at predicting electromagnetic decay amplitudes, we have to con- 
struct gauge-invariant expressions. Our Lagrangian involves massive scalar and 
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vector D-meson particles. We are faced with a serious complication, namely, 
the mixing of scalar and vector modes. It is a quite cumbersome enterprize to 
arrive at gauge-invariant expressions in the presence of such mixing phenom- 
ena [57|38ll39j . This is a known and non-trivial complication of the standard 
model where the Higgs boson may mix with the longitudinal component of 
the Z boson [ID]. A solution to this problem is to represent the 1~ D-mesons 
in terms of antisymmetric tensor fields [^42II43|H^ . The massive vector field 
is proportional to the divergence of the antisymmetric tensor. 

To proceed with this particular representation we demonstrate first that the 
results of [12], which were obtained using the conventional vector field repre- 
sentation, can be recovered with the tensor field representation. 

We start with the Lagrangian density, 

L = -{d,Dn (d-'D,^) + ^ D'^" Ml D,^ 
- ^ {{d^D,^) [$, D^^ - D,^ [$, {d.Dn} , (24) 

involving the kinetic term and its associated Weinberg- Tomozawa interaction. 
The antisymmetric triplet fields D^^ = —Dy^ and Z)^^ = Dj^^ with D^y = 
{D^^y, —D'^y, D^^^y) describe the heavy-quark multiplet partners of the fields 
D introduced in (JTj). Mi- denotes their mass matrix. Since the tensor field 
representation is not frequently used in the literature, we recall the definitions 
of the propagator and of the associated wave function. 



Ml 



(27r)4 F - Ml +ie 



X 



{Ml - kl 9ui3 + 9ixa K kf3 - g^p ky /c^ - (/i ^ v) 



(25) 



and 



P/.ev(p, A) -Pue^{p,\)\ , 
Eei(p,A)e.(p,A) = -^,. + ^, (26) 

A=l ^^^1- 

where the wave function is expressed most economically in terms of the conven- 
tional wave function e^(p. A) of a vector particle in the vector representation. 

To derive the on-shell scattering amplitude based on the interaction (!24l) . we 
have to reformulate with the tensor representation the technique developed in 



(0|D^,(0)|D(j>,A)) = e^,(p,A) = 
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pT] using the vector representation of spin one fields. Tlie on-sfiell part of the 
scattering amphtude in the vector representation takes the simple form, 



yon-shell ^ ^ ^iJP)^^^ yiJP)^-^ ^. ^) ^ (27) 
J,P 

where the projectors y^J'^q, q', w) are constructed to carry well-defined total 
angular momentum J and parity P. The projectors are polynomials in the 
initial and final 4- momenta of the Goldstone bosons, and g^, as well as in 
the total 4-momentum {w"^ = s). We are interested only in the = 
sector to generate the Z^*i(2460). We recall the appropriate projector 



2 V 

The merit of the projectors is their property of solving the Bethe-Salpeter 
coupled-channel equation analytically for quasi-local interactions. The partial- 
wave amplitudes M^-^^^s) are Lorentz invariant. They can be computed in 
terms of an effective interaction V^'^^^s) and loop functions J^'^^\s), 



jiJP)(s) = N^'^\s){l{s)-I{f,l,)}, (29) 

where the factors N^'^^'^ (s) reflect the presence of spin and angular momentum. 
The latter, if multiplied by an appropriate factor s", are polynomials in s and 
the masses of the intermediate states. The universal integral /(s) is defined 
in (fT5|) . The matching scale /im is taken to be the mass of the vector-meson 
ground state, /ia/ = 2012 MeV, following the same reasoning as given in 
Section 2.1. for the scalar case. 

We seek a set of tensor projectors, yjii^Jpiq, q', w) , with properties analogous 
to those of yli,^\q, q] w). They are defined in terms of the previous ones by 

yli^Jp {q,q;w) = \ y^^p {q.q]w)pp-\p^ y^^p (g, g; w) p„ 

-\p^l yy/^ (q^ Q; w)p(3 + j yiP (g, g; w) , (30) 
where p^ = Wfj, — q^ and p^i = — qfj,. By construction we have 



e^'^'iP) y%Z{q. ^) e'^^(p) = ^'e^''{p) y\f/\q. q; w) e^{p) , (31) 

where we made use of (1261) . The identity (13T]) provides the relation we are after 
for the scattering amplitude. 
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t;:-T' = E M^'''^ (^) yliZi^, g; w) , (32) 
j,p 

where the invariant partial-wave amphtudes are given by an equation of the 
form (l29il . The only modification compared to the expressions of [H] is a 
rescaling of the loop functions by a factor M^. The normalization factor reads 



N^^+\s) = ^M' + ^, ^s = ^M^+pl^ + ^m^+pl^. (33) 



Using a convention for the coupled-channel states analogous to f[T^ . the ef- 
fective interaction Vij^\s) implied by (l2ll) is obtained by a straightforward 
application of pT] , 



(1+) M^ + M^ (0+) , , ^ (M^ - M^) 2 2, ^ 

in terms of the matrix V^rp^^j{s) and the Cij coefficients defined in f|T7j) . The 
parameters M, M, m and m have the same meaning as in f|T6|) . 



In the limit M = M, we recover the expressions obtained in [12]. The invariant 
amplitude M^^~^\s) is identical to that of |T2] within a factor 2M^/3. This 
implies in particular that the predictions for the axial-vector spectrum are 
consistent with the expectation of the heavy-quark symmetry. 

The 1^*^(2460) state is generated dynamically from the interaction Lagrangian 
(l24l) . The scattering amplitude develops a pole at s = M^+. Close to the pole 
it has the form 



3 MM s - + i Ti+ MiH 



^rH^) ^ -T77TT7 . .r:.^;^ ^ (35) 



with the coupling constants (^j and the width parameter ri+. For the same 
reason as in flTSl) . we do not consider background terms. The normalization of 
the coupling constants, Qi, is such that in the heavy-quark limit they are iden- 
tical to those of ffTSl) . The values discussed in the following can be compared 
directly to those given in [T3] . 

To reproduce the empirical mass of 2459 MeV at leading order for the 1"*" 
state requires using an effective parameter /e// = 97.2 MeV. The coupling 
constants are then 



g^D' ^ 1.95, g^oD* ^ 0.049 
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gRODi ^ 2.25 , 



9K+D* ^ -2.25 



(36) 



Isospin-breaking effects in the coupling constants are again found negligible. 
For a quantitative study that considers chiral correction terms we refer to the 
end of this section. The values 0361) are identical (or very similar for g^ofj*) 
to those given in f[T^ . This approximate degeneracy is expected from heavy- 
quark symmetry. The strong width of the D*^{2460) is determined by Qt^od*- 
We get 55 keV using e = 0.010. Our result for the width parameter is five 
times larger than the value obtained in [15] based on the same interaction. 
The reason for such a discrepancy is again the neglect of isospin-breaking 
effects in the kaon and D*-meson masses in [15]. The uncertainties in the 
strong width discussed for the D*q(2317) apply equally for the D*]^(2460). 

2.3 Chiral correction terms 

It was shown in [13] that chiral correction terms at subleading order (Q^) pro- 
vide additional mass shifts for the D*q and D^^ states, such that their masses 
are consistent with the observed values when using / = 90 MeV rather than 
the effective decay parameters /e// of the previous sections. We construct 
these chiral corrections to the leading order interactions ([6]) and flM|) by ad- 
justing the expressions of [13] to the tensor representation of massive vector 
particles. There will be two types of contributions for s-wave scattering. On 
the one hand we include s- and u-channel exchanges of the D-meson ground 
states based on the leading order vertices involving a Goldstone boson and 
two D-mesons. On the other hand local 2-body counter terms (breaking chi- 
ral symmetry and chiral symmetric respectively) will be constructed. These 
different contributions are shown in Fig. [2l 




Fig. 2. Chiral correction terms from s- and u-channel exchange processes of the 
D-meson (left) and local 2-body counter terms (right). The thin and thick lines 
represent Goldstone and heavy-light mesons respectively. 

The s- and u-channel exchanges of the D-meson ground states at chiral order 
are generated by leading order (Q^) 3-point vertices, 
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C = i^ {d^, (9^$) (d'-D) - {d-D) {d^^) D^,] 

As derived in Appendix A, the decay of the charged D*-mesons imphes 

I^pI = 0.57 ±0.07. (38) 

The parameter in (1371) can not be extracted from empirical data directly. In 
the absence of an accurate evaluation within unquenched lattice QCD, the size 
of ^P can be estimated using the heavy-quark symmetry of QCD. As discussed 
in Appendix B, one expects at leading order 

9p = 9p- (39) 

The coupling constant gp contributes to the effective interaction V^^^\s) and 
via s- and u-channel exchange processes of the D-meson ground states. 
The specific processes involved are linked to the choice of representation for the 
massive vector particles. In the vector representation, the propagator contains 
a non-propagating O"*" component. In the tensor representation, it contains a 
non-propagating l"*" component. We derive the relevant contributions applying 
the on-shell reduction scheme of [TT] . 

In the scalar sector, only the u-channel exchange of the = 1~ contributes 
as a consequence of the tensor representation. We have 



2 Cu-ch j dx {q- q) 1^1- - {rfi^ - q ■ p) [m^ - p ■ q) 
q-q = ^Jrh^+pl,, ^fm^^Tj^ - pcmPcm x , 

q-p=-q-q + , p-q = -q-q+ , 

V~S=sjM^+ + + Plm = ^JM' + P'cn. + + Plm , (40) 

where (m, M) and (m, M) are the masses of the initial and final mesons. The 
parameter /ii- ~ 2059 MeV is taken to be the average of the vector D-meson 
masses. The coefficients C^^f^ are given in Table [1] for the channels relevant 
for the formation and decay of the 1^*0(2317). 
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Table 1 

The coefficients C^^'^^ that characterize the interaction of Goldstone bosons with 
heavy-meson fields for given isospin (I) and strangeness (S). The channel indices 1 
and 2 refer to the isospin basis used in |13j 



As already noted in [13] the influence of the u-channel process is of very minor 
importance for the formation of the D*q{2317). At / = 90 MeV and gp = we 
obtain a mass of 2304 MeV which is pulled down by 1 MeV only if we switch 
on gp = 0.57. 

The situation is different for the axial- vector state 1^*^(2460). Three processes 
contribute: the s-channel exchange of the 1~ state and the u-channel exchanges 
of the 0~ and 1~ charmed mesons. The two u-channel contributions give 



1 



's) = / ^ {a, (1 + x') + p,^ X (1 - x') A,} 



u—ch \ 



f 



+ ((m - 3 nl-) + (/i?- (2 s - m) - u^) 

+ (M2 _ ^2 _ ^) ^ ^2_ ) ^2 ^ (^2 ^ 

+ (M^ + m) ((M^ -m) + m(-M2 + 2s + m))} , 

u = M"^ + M'^ - s + 2q - (41) 

where the kinematics is given by (140|) . We use /iq- = 1918 MeV as the average 
mass of the 0~ charmed mesons. The influence of the u-channel exchange 
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interaction (HT]) on the formation of the D*2(2460) state is somewhat more 
important than it is in the scalar sector. For / = 90 MeV and Qp = Qp = 0, 
we obtain 2441 MeV, a mass which is pushed up by 5 MeV upon incorporation 
of (HTi) with gp = gp = 0.57. The effect is dominated largely by the exchange 
of the 1~ state. 



We turn to the s-channel exchange. In contrast to [13] based on the vector- 
field representation of the spin-one D-mesons, there is a contribution from the 
s-channel exchange of the 1~ state within the tensor-field approach, 



Vn (^) = - + - + ' (42) 

where the coefficients Cfi'f^ are given in Table [H The combined effect of 
dUD and (02]) yields a resonance mass of 2432 MeV for / = 90 MeV and 
gp = gp = 0.57. 

We turn to the local 2-body interaction terms, considering successively the 
chiral symmetry breaking and the chiral symmetric terms introduced in [T3] . 

At chiral order Q^, the following terms break chiral symmetry explicitly. 



C^.SB = -2ciDxoD-{4co-2 ci) (D D) tr xo 

+ ^^^^5tr(<|.xo$) + ^/^{$,{$,Xo}}5 

■ ci D^p xo + (2 5o - 5i) 5"^) tr xo 

2 Co - Cl ^ -^g ^ ^\ Cl 



Z}«^5"'3tr(<|.xo$)-^/^./3{<^>,{<^>,Xo}}5"^ (43) 

where the matrix xo is defined in ([8]) . The parameters c\ and c\ are determined 
by the empirical mass differences of the = 0~ and = 1~ charmed 
mesons. According to [13] we have 



ci~0.44, gi~ 0.47. (44) 

The parameters Cq and Cq could in principle be determined by unquenched 
lattice QCD simulation upon studying the pion- and kaon-mass dependence 
of the D-meson ground states. So far they are unknown. We construct the 
effective interaction 
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C j > C j 



+ ■••, (45) 



where the dots represent higher order terms that we neglect. The coefficients 
Cl^fi^^ , Cl^^f^ , Cj^'P and C^'P are given in Table [H 

The chiral symmetric terms of order-Qy contributing to s-wave scattering are 



CcT = ^^y^ DDtT ((9,,$) id''<^)) -^D id^<^) D 



^ + D^^ D^^ tr ((9^$) (9''<l>)) + D„;3 (9^$) (9^$) .(46) 



They imply the effective s-wave interactions 



>(s) = VXl 'is) + ■ ■ • , (47) 

3MM ^ 

where the dots represent again higher order terms. The coefficients C2 ' and 
C^^'^^ are given in Table [TJ The full effective interactions V'^^'^^s) and V'^^'^^s) 
are then defined as the sums 

V{s) = VwTis) + Vu-ch{s) + Vs-chis) + V^.sb{s) + Vcris). (48) 

The number of unknown parameters, Co,2,3 and £0^2,3 appears large at ffist. A 
free fit to the masses of the D*q{2317) and D*^(2460) states only would not 
be significant. Additional constraints from QCD should be used. According to 
|13] . the parameters q and Cj are degenerate in the heavy-quark mass limit, 
i.e. we expect 



Ci ~ Q . (49) 

It is reassuring that the values for ci and ci given in fl44p are quite compatible 
with the expectation of the heavy-quark symmetry relations f H9|) . We consider 
further constraints from QCD as they arise in the limit of a large number of 
colors Nc [l5]. Since at leading order in a 1/iVc expansion single- flavour trace 
interactions are dominant, the corresponding couplings should go to zero in 
the Nc 00 limit, suggesting 



Cl C3 _ Cl „ C3 

Co ^ ^ , C2 ~ -y , ^0 - Y ' ^2 ^ ~ 2" • 



20 





0+ 
(0,1) 


0+ 

(0.1) 


(io) 


1 + 

(0:1) 


1 + 

(0,1) 


1 + 

(io) 


Mr, \Mo\l] 
iVi^ [lVi6 V J 


9'^1 7 fi 
Zol 1 .0 


9'^1 7 fi 
Zol 1 .0 






941^0 9 


ZOOo 


r fMeVl 

i- 1 iViC- V 


n 14 


95 


2.18 


14 


25 


1 8 


loi 1 

li/i 1 


3.27 


3.27 


0.24 


2.97 


2.97 


<0.05 


152 1 


2.50 


2.50 


1.37 


2.42 


2.42 


1.4 








2.12 






2.5 


e 


0.01 


0.02 





0.01 


0.02 






Table 2 

Masses, widths and coupling constants for dynamically generated O"*" and 1"*" states. 
The coupling constants gt corresponds to specific decay channels labeled in the 
isospin basis used in [13]. We use / = 90 MeV, gp = gp = 0.57, C3 = 1.0, C3 = 1.4. 
The remaining parameters are implied by the relations (1441 [49l [50]) . 

In the combined heavy-quark and large- A^^^ limit fl^ ISUl) . we are left with 
one free parameter only, C3 = 63. The optimal value 03 = 03 = 1.2 (together 
with the parameters ci and ci given by (l44l) . gp = gp = 0.57 and / = 90 
MeV) predicts 2330 MeV and 2448 MeV for the masses of the D;o(2317) and 
D*]^(2460) states respectively. 

For the determination of the width parameters, it is important to reproduce 
the masses accurately. We allow therefore for small variations of the parame- 
ters around the heavy-quark scenario, leaving the large- A^c relations ( l50l) un- 
touched. A precise reproduction of the scalar and axial- vector state masses is 
achieved with C3 = 1.0 and C3 = 1.4. Detailed results are collected in Table [2l 
For the mixing angle e = 0.01, the strong decay widths of the D*q(2317) and 
Z^*]^(2460) are both 0.14 MeV. These widths are significantly larger than the 
leading order estimates discussed above. This effect is mainly a consequence 
of somewhat larger coupling constants g^^jj^ = 2.50 and g^^j^} = 2.42 given in 
Table [2] (where they are denoted by g2). The sensitivity to the mixing angle 
e is illustrated by the 3rd and 6th row. Taking e = 0.02 instead of e = 0.01 
leads to a strong width of 0.25 MeV rather than 0.14 MeV. 

We briefly comment on the previous results of [13]. In that work a different 
scenario was investigated. Using the conventional vector-field representation 
of the 1~ charmed mesons, it was assumed that the axial-vector resonance 
Dl{24:20) was a member of the exotic sextet, predicted at leading order by 
chiral coupled-channel dynamics [T^IITB] . The crucial question is whether chi- 
ral correction terms reduce the weak attraction predicted at leading order or 
possibly enhance it. This is clearly sensitive to the implementation of chiral 
symmetry. For example, in [T6IT7] based on the broken SU(4) flavour symme- 
try, the sextet states are weakly bound only and therefore quite broad. In [13] 



21 



chiral correction terms were tuned in such a way as to pull down the exotic 
axial state with (I,S) = (|,0) to match the properties of the D*(2420). The 
invariant ttD and tiD* mass distributions as measured by the BELLE collabo- 
ration [in] were used as an additional constraint. It was argued that the scalar 
heavy-quark partner of the exotic axial state decouples from the ttD channel 
and therefore is not seen in the data. Based on the large- iVc relations, that 
were not considered in [13], we would deem this scenario unlikely. In order to 
compare to the large- iVc scenario advocated in the present work we included 
in Table [2] the rows 4 and 7, which give the characteristics of the exotic (|, 0) 
states. Like in [13] the scalar state is quite narrow with a mass below the rjD 
threshold. We obtain a mass and a width of 2411 MeV and 2.2 MeV respec- 
tively. As seen from the coupling constants in Table [21 the state couples most 
strongly to the rjD and K Ds channels. 

In Fig. [3] we confront the imaginary part of the nD nD amplitude with 
the invariant tiD mass distribution measured by the BELLE collaboration 
[l6] . Such a comparison is approximate since it does not resolve the structure 
of the initial state. The empirical distribution is dominated by the broad 
(|,0) state, a member of an antitriplet like the D*q(2317), and the tensor 
state 1)2(2460). The contribution of this state is illustrated by the histograms. 
The possible presence of a narrow (|,0) state is not excluded by the present 
data. It is interesting to observe that the exotic state leads to a dip in the 
mass distribution rather than a peak. This is a consequence of the nearby "qD 
channel that couples strongly to that state. With the exception of a strong 
cusp effect at the KD threshold in the (0, —1) sector, there is no further strong 
signal of any sextet state in this scenario. 

A striking prediction of the large-A'^c scenario is a clear measurable signal 
of an exotic axial state in the rjD* invariant mass distribution with a mass 
of 2568 MeV. It lies above the rjD* threshold and has a width of about 18 
MeV. In Fig. [3] we confront the imaginary part of the ttD* ttD* amplitude 
with the invariant irD* mass distribution measured by the BELLE collabo- 
ration [l6]. The empirical distributions show two axial and one tensor states. 
The contribution of the tensor state is again illustrated by the histograms. In 
conventional approaches the tensor state 1^2(2460) is grouped together with 
the Z)j|^(2420) state to form a heavy-quark multiplet. Within the hadrogen- 
esis conjecture we would expect to generate that multiplet dynamically via 
coupled-channel effects once the light vector mesons are considered as addi- 
tional and explicit degrees of freedom. The theoretical amplitude of the present 
work describes only the broad state, which has a width of about 300 MeV. In 
contrast to the ttD distribution shown in the left panel of Fig. [31 we do not 
predict any significant signal in the ttD* distribution that one may use to dis- 
cover the exotic axial state. This reflects a coupling constant of that state to 
the irD* channel that is almost compatible with zero. Nonetheless, the exotic 
axial state could be discovered by ongoing experiments once the invariant r]D* 
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Fig. 3. Mass spectra of the 0)-resonances as seen in the vr L'(1867)-channel (l.h. 
panel = 0"*") and vr D*(2008)-channel (r.h. panel = 1"*^). The solid lines show 
the theoretical mass distributions. The data are obtained from the B 7rD(1867) 
and B ttD*{2008) decays |46j . The histograms indicate the contribution from 
the J = 2 resonances 1)2(2460) as given in |46j . 

mass distribution is analyzed. The discovery of the scalar state would require 
a measurement of the nD invariant mass with an energy resolution of a few 
MeV as may be possible with the PANDA experiment at FAIR. 

We close this section with the remark that very similar results are obtained 
with the vector-field representation once the parameters are chosen in ac- 
cordance with the expectation of large- A'c QCD rather than with the fitting 
procedure recalled above p^flTj . 



3 Electromagnetic interactions 

To compute the radiative decay of the scalar and axial-vector states gene- 
rated in the last section we need to couple the photon to the hadronic fields, 
ensuring consistency with the U(l) gauge symmetry of electromagnetic inter- 
actions. We consider first the electromagnetic interactions obtained by gauging 
the hadronic interactions of Section 2. We add terms involving Goldstone and 
D-meson fields which are separately gauge-invariant and of chiral order Q^. 
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Such terms are not sufficient to reproduce the available data satisfactorily. We 
introduce an additional 4-point vertex of order to evaluate the need to in- 
clude higher chiral orders. With such a term, the data on the radiative decays 
of the D*q(2317) and D*^(2460) mesons can be understood at the expense of 
an unnaturally large value for the coupling constant of the term of order Q^. 
We expect this result to reflect the effectiveness of the underlying theory and 
more precisely the absence of light vector mesons in the dynamics of radia- 
tive transitions. We introduce them as additional explicit degrees of freedom 
and will indeed find that this leads to natural values of all the parameters, 
justifying a posteriori this extension of the theory. Light vector mesons are 
not included in the coupled- channel part of this work but should become part 
of that scheme in a later stage to study the formation of tensor states, in 
particular the low-lying 2^ state expected at 2573 MeV. 

We gauge the hadronic Lagrangians (El [241 EZD using the covariant derivatives 



d^^ + ie[Q,^]A^, d^D + ieDQ' A^, (51) 

where e = |e| and the charge matrices are defined by 



Q 



^ 

0-1 







'^O 0^ 



Q' 



3/ 



1 
1 



(52) 



The radiative decays of the D*q(2317) and D*^(2460) mesons will involve only 
3-point and 4-point vertices. Gauging 4-point vertices, such as the Weinberg- 
Tomozawa terms in ([6]) and flMl) or the chiral correction terms P6l) , would lead 
to 5-point vertices of no relevance. The kinetic terms in ([6]) and flM]) imply 
the 3-point couplings 

C,.^.=i'^Ani ((9^$) [g, + ieA^^[DQ\d^D)-{d^D)Q' D) 
-leA^ D^, Q' (d^Dp,) + ieA^ {d^D'^') Q' D^, . (53) 

Given the interactions ([HI [Ml [351) only, the electromagnetic decay amplitudes 
0"'' — > 7 1~ and 1+ ^ 7 0~ or 7 1~ are zero identically as already noted in |48j . 
These decay processes probe on the other hand the 3-point hadronic vertices 
introduced in (pTj) . The importance of these terms in determining decays con- 
trasts with the minor role they play for the formation of the D*q(2317) and 
Z>*j^(2460) states. The gauging of the interaction (1371) yields the terms 
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+ ^ {d,, id''^) A''Q'D + A'D Q' (9^$) D^,} 

+ ^ ^ e'^""^ {/^M- [Q, A. (d^D^ts) + (S^/^.^j) [Q, $ ] A^ D^,} 

- t e^'^"^ {d^, A- g'D.^ - A-^ D^p Q' D^,} . (54) 

These terms are part of the gauge-invariant vertex associated with a leading 
chiral power Qx? hence showing that the standard counting scheme requires 
the photon field to be of order Q^. 

We construct the chiral correction terms of order Q"^. They are gauge- invariant 
separately as they involve the electromagnetic field strength tensor F^,^ (of 
chiral order Q"^)- Terms of a similar structure were discussed previously, for 
example in [IS] , using a Lagrangian where massive spin 1 fields were described 
in the vector representation. Clearly the Lorentz structure of these terms will 
be different in our case. 

Chiral correction terms of order describe anomalous processes. A photon 
hitting a pseudoscalar charmed meson may convert the latter into its heavy- 
quark partner, a charmed vector meson. In the absence of such processes the 
decay amplitudes — > 7 D and D* ^ Dg would vanish identically, contrary 
to experiment. The leading order anomalous vertex should be of the form 

A.m. = ^ e^^^p {{drD^^) {ec + eg Q) (d^D) 

+ {d''D)[ec + eQQ){drD^^)}. (55) 

The vertices f l55l) carry the leading chiral power [^49|50j . The radiative 
decay properties of the charmed vector mesons suggest (see Appendix A) 



CQ = 0.91 ±0.10, ec = 0.13 ±0.05. (56) 

The values of eg and ec given in ( l56l) reproduce the empirical branching ratios 
of the D*~ ^ L>~ 7 and D*^ j decays. The data on the D* ^ D^'j 

and D* DsH^ suggest the smaller value cq ~ 0.52 (see Appendix A). The 
anomalous interaction (!55l) by itself is at odds with the heavy-quark symmetry 
of QCD which relates the interactions of pseudoscalar and vector D-mesons. 
Additional terms are required for consistency to parameterize the magnetic 
moments of the charmed vector mesons. They read 
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A.m. = ^ F^^'' (d'^D^^) (e Q' - ec + eg Q) {d^Dp,) . (57) 

The parameters ec and eg are related to the magnetic moments of the charmed 
vector mesons by [^42|I43] 



M^j . 2 ^ _ Mz,; . 1 



2 V - 3^7' 2 
Md* / 1 \ 



At present there is no empirical information on these magnetic moments of 
help to determine the values of ec and eq. Note that the term proportional to 
Q' in ( 1571) cancels a corresponding contribution which is implied by minimally 
gauging the kinetic term ( l53l) . 



It is instructive to interpret the result ( l58i) in terms of the constituent quark 
model. The contribution from ec reflects the magnetic moment of the charm 
quark. It is SU(3) flavour-blind. In the heavy-quark mass limit, the parameter 
e-c approaches a constant. In contrast, the term proportional to eg models 
the contribution of the magnetic moment of the light quark. It is flavour- 
dependent, being proportional to the charge matrix of the light quarks Q. In 
the heavy-quark mass limit, eq scales linearly with the charm quark mass. 
To leading order, the parameters ec and eq are related to ec and eq by (see 
Appendix B) 

= cq , ec = ec. (59) 

We emphasize that the two sets of terms ( l55i) and ( 1571) are the only chiral 
corrections to the leading order Lagrangian at order Q"^. It does not appear 
possible to fit the data on the 1^*0(2317) and D*^{2AQQ) radiative decays with 
these electromagnetic terms together with the Lagrangians ( !53|) and (15^ (see 
Section 6 for details). This suggests that higher order terms should be consi- 
dered. We introduce an additional term involving one Goldstone boson field 
and proportional to F^^, and to the charge matrix Q of the light quarks (ex- 
pected to be dominant in the heavy-quark mass limit 



+ (9"D,;3)[(<9"$),g] D^r)]. (60) 
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The heavy-quark symmetry predicts the relation (see Appendix B) 

ep = ep . (61) 

One expects \ep\ ~ e ~ 0.303. According to standard counting rules the 
vertices (1601) carry the leading chiral power Q^. It should be emphasized, 
however, that the counting rules depend on a naturalness assumption, i.e. 
that the dimensional parameters of the effective Lagrangian density scale with 
appropriate power of the chiral symmetry breaking scale 47r/ ~ 1131 MeV (or 
more pragmatically with the mass my of the lightest degree of freedom that 
is integrated out). This is the rationale behind the particular representation 
of the vertex (!60|) in terms of the dimensionless parameters ep and ep. One 
may assign the vertex (160|) the order Q^. On a formal level that would justify 
to consider the effect of (1601) while neglecting additional hadronic vertices of 
chiral order Q^, such as the SU(3) flavour breaking effects in the coupling of 
the Goldstone bosons to the D mesons (see fl37j) ). 

There is no empirical estimate available of the size of ep and ep, in particular 
no data on the three-body decay process D*~^ — > 7 tt"'" sensitive to ep. These 
parameters have therefore to be determined indirectly. We shall see in the 
section on numerical results that the magnitude of the parameter ep is larger 
than expected from the naive naturalness assumption if we try to reproduce 
the data with the effective Lagrangian derived above. Such a phenomenon is 
not unusual in effective field theories, though asking for a physical explanation. 
We assume that the size of ep is a manifestation of the need to incorporate 
additional degrees of freedom, the most natural ones being the light vector 
mesons. Once they are incorporated one expects the fitted value of ep to 
reduce significantly. We take this as a motivation to explore the role played by 
the light vector mesons in the radiative decays of the Z)*q(2317) and D*-^^{24Q0) 
states. 



3.1 Light vector mesons 



For the light vector mesons, there is no well-defined chiral power counting 
scheme to the authors' knowledge, at least for the case of virtual vector mesons. 
The treatment of these degrees of freedom is therefore more phenomenological 
than the chiral expansion approach used for Goldstone bosons and heavy-light 
mesons. The 3-point electromagnetic vertex probed when considering the light 
vector mesons is the anomalous process analogous to ( l55i) where a photon 
excites a Goldstone boson into a vector meson. It reads 

>Ce.m. = e^^^p tr((a"$) [Q, ^rV^^]^ , (62) 
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where we use the representation 



V, 



UJ 



V 



V2K 



UJ 



(63) 



The particular form (162|) is equivalent to the anomalous photon coupling used 
for example in [51] . From the radiative decays of the light vector mesons con- 
sidered in Appendix A, one derives conflicting values for the parameter e^. 
A quantitative description requires the inclusion of SU(3) flavour breaking 
effects. Only three processes, K*^ — > K^'^, K*^ — >■ K^'^ and — > 777, are rele- 
vant for the radiative decays of the D^q and D*^ states. We will include SU(3) 
breaking by using the phenomenological couplings obtained in Appendix A 
for each of these processes, 



|^(i^*-.i^o7)| ^ 0.119 ±0.006 



I ( 777 ), 

\^A I 



0.053 ±0.001. 



Va I 



0.090 ±0.004. 



(64) 



The anomaly vertex ( l62l) can contribute to the radiative decays of the D*- 
mesons only in the presence of additional hadronic 3-point vertices. The light 
vector meson that is created by the photon must be absorbed by a heavy 
meson. The corresponding interaction is constructed by analogy with ([3j 



9y_ 
2/ 



9v 



4/ 



(65) 



The electromagnetic terms implied by the covariant derivatives f l51l) should 
be added. The particular form of (165!) was guided by the Lagrangian (137|) and 
the flavour SU(4) limit. We complement (!65|) by additional electromagnetic 
interactions analogous to (!60l) . 



ey 
2/m2 



V 
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^ / ''V 

- 7^ (daD) {V^,, Q} (d'^D) 

+ 4^ F^- (d^D^n {V, Q} {d^Dpr) . (66) 
/ ^y 

The hadronic vertices ( 1371 1651) should also have some impact on the forma- 
tion of the molecules based on the leading orders coupled-channel interaction. 
The influence of the parameters gp and gp was studied in [13] and in Section 
2.3. Their effects on the mass were found to be rather minor, confirming the 
expectation that the corresponding interactions are of subleading order in a 
chiral expansion. Nevertheless, this issue deserves further studies. In particu- 
lar, the role played by additional inelastic channels involving the light vector 
mesons should be worked out. We anticipate that they are responsible for the 
formation of tensor molecules. We will return to this issue in the next section. 

The various parameters in ( |65l [66l) are correlated by the heavy-quark symme- 
try of QCD. As shown in Appendix B, one expects 

9v = gv, 9t = gr, 

ey = ey , ct = ct, ge = , (67) 

at leading order. The size of gv and gv can be estimated by the phenomeno- 
logical assumption of universally coupled light vector mesons and the KSFR 
relation to be 



with the universal vector coupling constant \g\ ^6 [52j. The phase of the 
coupling constant gv is not determined by this assumption. 

The estimate of the size of the coupling constants gx and is more difficult. 
While lacking QCD lattice simulations, we may use a flavour SU(4) ansatz, ad- 
mittedly a rather rough and questionable tool. In order to do so, we introduce 
flavour SU(4) multiplet fields 
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^[16] 



[16] 




(69) 



in terms of the SU(3) multiplet fields $, D and V^,y, D^^. We construct a 
minimal SU(4) invariant Lagrangian density that comprises the interaction 
terms introduced in fl371 [65D. 



5(7(4) 



+ ^tr (9.y[r6l)e'l(5/.^?e,)-^ 
The ansatz (ITOl) suggests the identification 



tr(9^$[i6l)Vf,';j (9.$[i6])(.70) 



9t = 9e = 9p = 9i 



9v = 9t = 92 



9v = 9p = 93- (71) 



The result fl7ip deserves and requires some discussion. We observe that ac- 
cording to fl39l [67j) the coupling constants gi with z = 1,2,3 in flTOj) should 
approach a universal value in the heavy-quark mass limit, for instance the 
parameters gp and gv- Confronted with the values fl38|) and fl68l) . this expec- 
tation appears verified within less than 20%. However, the result flTTl) would 
predict also a common value for gp and gx- We point out that this implication 
is troublesome: whereas gp approaches a finite value, the parameter gx must 
vanish in the heavy-quark mass limit. This is an immediate consequence of 
the fact that the QCD action is linear in the charm quark mass. While the 
term proportional to gp involves one derivative acting on a heavy field, the 
term proportional to gx involves two. It is therefore clear that we must not use 
the relations (ITTil as they stand. From (TrTI) we retain only the predicted phase 
relations for the coupling constants, i.e. we assume all coupling constants to 
take positive values. In addition, for the unknown parameter gE, we would 
anticipate the range 



9p <9e < 9v 



(72) 



a conjecture consistent with the expected scaling behaviour at large charm- 
quark masses. The estimate of the remaining parameters gr — gr is most 
uncertain. From 0711) we would expect the range 



my 

< Qrp < —— Op , 



(73) 
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with my = 776 MeV and My = 2000 MeV representing the typical masses of 
hght and heavy vector meson. 



3.2 Radiative decay of molecules 

We display in Fig. H] the various graphs contributing to the radiative decay 
of the 1^*0(2317) and D*^(2460) states which are linear in the hadronic three- 
point vertices. The D^q and mesons generated by coupled-channel dyna- 
mics appear as blocks, the sum of the bubble diagrams defining the effective 
propagator of the resonance state. It is drawn as a double line. The dynamics 
of the decay is contained in a single loop. Solid lines represent the propagation 
of the pseudoscalar or vector mesons. The thin lines stand for the light mesons 
and the thick lines for the heavy mesons. The wavy line is the photon. 

We display in Fig. [5] two-loop diagrams which, at first sight, should be relevant 
for the radiative decay width of scalar and axial- vector molecules. The figure 
shows contributions where the photon couples to the resonance directly. We do 
not resolve the structure of the vertex since these diagrams vanish identically. 
This is an immediate consequence of using the tensor-field representation of 
the spin one particles. Consider for instance the decay of a scalar molecule, 
where the final state carries = 1~ quantum numbers. All diagrams in 
Fig. factorize into two contributions. Since the left part is contracted with 
the antisymmetric wave function of the 1~ state, it vanishes identically. This 
contribution carries indeed the two indices a and (3 and depends on one 4- 
momentum only, the 4-momentum of the final state. It must be symmetric in 
the indices a and j3 and therefore vanishes if contracted with the wave func- 
tion. Analogous arguments hold for the decay modes of axial-vector molecules. 
There are two possibilities. If the final state carries = 0~ quantum num- 
bers, the antisymmetric and outgoing indices of the intermediate axial- vector 



Fig. 4. Single-loop diagrams contributing to the decay amplitude of a scalar or 
axial- vector molecule (see text). 
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Fig. 5. Two-loop diagrams contributing to the decay width of scalar and axial- vector 
molecules where the photon couples to the resonance directly (see text). These terms 
vanish identically in the tensor representation of massive vector particles. 

propagator are contracted necessarily with the 4-momentum of the final state. 
If the final state carries = 1~ quantum numbers, the final 3-point vertex 
involves the antisymmetric tensor e^i^ap due to parity conservation. At least 
two indices must be saturated by the 4-momentum of the 1~ D^-meson. 

A striking conclusion of our discussion is that the radiative decay of a scalar 
or axial- vector molecule is determined fully by an effective hadronic resonance 
vertex. The detailed structure of the resonance propagator is not probed. We 
emphasize that this is a consequence of using the tensor-field representation 
of spin-one particles. The application of the vector-field representation would 
require a detailed study of the resonance propagator. This could be quite 
cumbersome. In particular, enforcing gauge invariance is highly non-trivial. 
Our observation implies that even in the presence of the additional coupled- 
channel interactions fl37l [651) . the formal evaluation of the diagrams of Fig. H] 
would be unchanged. 

Because of the particular structure of the decay diagrams, it is useful to in- 
troduce flavour triplet fields, R and R^^, which interpolate the scalar and 
axial-vector molecules. All what is needed are the coupling strengths of the 
molecules to the hadronic final states. We introduce the effective coupling 
constants and g^, 

^i?, = -0 {id,D) {d'^<^) R + R {d^<^) {d,D)} 

- ^ {D'"' (d,^) {d^Rru) - {d^Rru) {d^<^) D^""} , (74) 

where we do not write explicitly terms linear in the photon field that are re- 
quired by gauge invariance. We may add terms involving the electromagnetic 
field strength tensor F^j^ that are gauge-invariant separately. Like the parame- 
ters gR and gR, the structure and size of these terms have to be extracted from a 
coupled-channel computation. As will become clear in the subsequent section, 
the role of such terms is very minor. Within our renormalization scheme such 
contributions vanish identically for the particular choice where the matching 
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scale jj^M and the mass of the hadronic final state are degenerate. This de- 
generacy is almost realized. For on-shell conditions, the first vertex of (!74|) is 
equivalent to a vertex of the generic form D ^ R as implied by the projector 
technique described in Section 2.1. This follows from the replacement 



2 {d,,D) (S'^^) R^D<^ {d^df,R) - {d^df,D) <i>R-D (9^9^$) R . (75) 

In general the two vertices would give different results in one-loop diagrams. 
Within our renormalization scheme discussed in the next subsection, the ver- 
tices are equivalent. The corresponding decay amplitudes coincide up to terms 
that may be generated by effective molecule vertices involving F^j,. As ar- 
gued above such terms vanish identically for the particular choice where the 
matching scale fiM is identified with the mass of the hadronic final state. An 
analogous argument shows the equivalence of the second vertex of flTll) with 
the generic form (5^-0"^) $ {d^ Rpr) implied in Section 2.2. This vertex leads 
to the form (l27|l of the scattering amplitude for = l"*". 

At leading order the magnitude of the parameters gn and qr are determined by 
the coupling constants as extracted from the pole structure of the scattering 
amplitude (see f|T8 l [35]) ). We identify 



^ 4V2/M0+ (0+) ^ 4^/6/Mo+ (0+) 
MoV -Ml -ml MoV - Ml^ - ^''^^ ' 

„ _ 4V2/Mz)> (1+) _ 4V6/M^; (1+) 

" Ml^ - Ml, - m\ " Ml^ - Ml. - ml ^^^^ " ^ 

As detailed in Appendix B, heavy-quark symmetry predicts 



gR = 9r, (77) 

at leading order. While fl77|) is realized quite accurately, the SU(3) fiavour 
breaking effects in the coupling constant gn and gf> are sizeable as can be 
inferred from the values given in Table O The relation (177|) is satisfied at the 
10 % level 



— ~ 986 ^^^^ ~ 989 (78) 
In contrast, the SU(3) relations 

gn ^ 0.720 ^ 1-710 g^^ , ~gR ^ 0.710 g^^^l ^ 1.691 gi^^} , (79) 
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are violated at the 95% level. The value qr ^ 4.28 derived from the r]Ds 
channel with / = 90 MeV is almost a factor of two larger than the value 
obtained from the KD channel. 

Before turning to the renormalization issue, we discuss additional resonance 
vertices involving light vector mesons. The latter arise necessarily in a coupled- 
channel computation once any of the parameters gviQVigTiQT oi Qe is non- 
vanishing. Though we are not presenting results of such a computation, we 
study the relevance of light vector meson in radiative decay processes of scalar 
and axial-vector molecules. We introduce the effective resonance vertices (lea- 
ving out additional terms linear in the photon field required by gauge invari- 
ance) 



where the parameters gn-, Qh and qh are unknown at this stage and degenerate 
in the heavy-quark limit. As worked out in Appendix B, we have indeed 



For the axial-vector molecules, additional vertices that are not on-shell equi- 
valent to those of flHOl) may be constructed. Such terms are suppressed by 
their d-wave phase-space behaviour. The presence of the vertices (1801) leads to 
additional diagrams contributing to the radiative decay amplitudes of scalar 
and axial- vector molecules. Such diagrams are part of Fig. Hlgiven the meaning 
of the solid lines which represent the propagation of either pseudoscalar or 
vector mesons. 

3.3 Renormalization 

We finally turn to the renormalization issue. The decay diagrams are ul- 
traviolet divergent. Applying the Passareno-Veltman reduction [53], which is 
rigourously justified within dimensional regularization, the integrals of Fig. H] 
may be expressed in terms of a set of scalar integrals of the form 




^ [r [d^Vru) [d^Dn + {d,Dn [d^Vr.) R} 
jj {{d^R,u) (drVn D-D (drVn {d'R,u)} 
f- e-^^^ [{dPR,,) (a^K.) + D^p (9-V;.) [d^'K,]], (80) 



Qh = 9h = 9h ■ 



(81) 





1 
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r d'^l - r dH 

hb = J Sail) S^il + p) , lab = J Sail) Sf^il + p) , 

Jabc = +^ / ^ ^'^(0 S,il + q) S^il + p) , 

Jabc = +^ / S^i^) Sbil + P) S^il + P) , (82) 

where we focus on the physical hmit with space-time dimension four. In our 
convention the 4-momentum p^ characterizes the decaying molecule, is the 
4-momentum of the photon and Pfj, = Pfj, — Qfj, is the 4-momentum of the 
hadronic final state. The quantities defined in (11591) are evaluated in Appendix 
C. 

Divergent structures arise only from the tadpole integral and the two- 
propagator integrals lab and lab- Since the combinations 



lab — lab , lab 5 -q ! (83) 

— 

are finite as space-time dimension approaches four, one may take the viewpoint 
that all divergent structures are caused by the tadpole integral la (see also 
[54]). 

Any divergent contribution can not be renormalized by simply adding a counter 
vertex where the photon couples directly to the resonance. Since the resonance 
is formed by coupled-channel dynamics its radiative decays amplitude must 
be renormalized by a loop subtraction mechanism similar to the one intro- 
duced in This is illustrated in Fig. O Typically the counter terms needed 
are appropriate 5-point vertices as indicated in Fig. [6] by an solid square. A 
contribution to the radiative decay amplitude is generated by contracting two 
identical lines giving rise to a tadpole-type contribution. The latter have to 
cancel the tadpole contributions from Fig. H] as discussed above. 

There is yet another important observation to make. In Fig. H] the integral 
arising from the second and fourth diagrams have a structure similar to the 

Fig. 6. Examples of counter loop contributions that renormalize the decay amplitude 
of molecules. 
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loop functions that build up the resonance state. Therefore the renormaliza- 
tion of the coupled-channel dynamics and of the radiative-decay amplitude 
are necessarily related. Since the coupled-channel dynamics involves an infi- 
nite number of Feynman diagrams, one has to leave the well-trotted path of 
perturbative renormalization, i.e. infinite sets of counter terms present in ef- 
fective field theories have to be introduced at each order. We emphasize that 
one should carefully discriminate two issues. First, the renormalization scale 
independence of a scattering amplitude and second, a possible scheme depen- 
dence of how the infinite number of diagrams are treated. We point out that 
the application of the on-shell reduction formalism developed in [^33|llj . 
as is implied by the unique existence of an algebra of covariant projectors, 
may be considered as a scheme dependence. The additional ingredient, the 
requirement of a smooth matching of scattering amplitudes unitarized in dif- 
ferent channels, can be considered again as a scheme dependence. It is an 
economical, though not unique, procedure to build crossing symmetry into 
the scheme. From this point of view the matching parameter fiM in f|T^ [291) 
reflects a scheme dependence, not a renormalization scale dependence. 

As a consequence of the on-shell reduction formalism our coupled-channel 
dynamics is based on, neither the effective potential V^'^^\s) nor the loop 
functions J^'^^\s) gain any contribution at leading order from a tadpole-type 
loop integral la- Using the leading order two-body interaction as the driving 
force of a Bethe-Salpeter equation would introduce plenty of tadpole contri- 
butions [53PT] . As was discussed in great detail in PSIITT] . the contribution of 
reduced tadpoles can be trusted only at a level where one computes one-loop 
corrections to the Bethe-Salpeter interaction kernel, i.e. for two-particle irre- 
ducible diagrams. While for the latter conventional power-counting arguments 
are applicable, they are not for the reducible contributions, i.e. those summed 
by the unitarization. 

Closing our chain of arguments, we arrive at the result that reduced tadpole 
contributions in the decay amplitude should be dropped in our leading order 
computation. In addition a finite renormalization is applied to the integral lab- 
It is identified with the expression introduced in f|T5l) with 



hb = I{p^) - /(/^L) > ma = M , mb = m, (84) 



and the matching scale hm used in the coupled-channel computation. 
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4 Radiative decay of scalar molecules: D*q{2317) ^ D*{2112) 



We express the partial decay width ro+-^^i- of the D*Q(2317)-meson in terms 
of the transition amplitudes M^f^^-^_ . Due to gauge invariance, the transition 
tensor is characterized by a single number do+^-yi- defined by 



= 4 (g, A,) eif^ip, A,) {g^^ {p ■ q) - g"} ^ , (85) 

where e^(g, Xg) and eapip, Ap) denote the wave functions of the outgoing photon 
and vector meson and p^ = Pi^ + Qfi the decaying resonance momentum. The 
radiative width reads 



47r I 2 Mn+ 



^6) 



with = p"^ and Mq+ = p^. The decay parameter can be obtained by the 
projection formula, 



since the decay amplitude is transverse with respect to the photon momen- 
tum and antisymmetric in the indices a ^ (3. This implies that the off-shell 
amplitude ppM""^'^ is orthogonal to and pa-, i.e. is characterized by two 
parameters only. With (!87|) we project onto the relevant component. 

There are 5 classes of contributions as shown in Fig. HI All terms have the 
topology of a one-loop Feynman diagram. A given contribution is either pro- 
portional to the coupling constant or gu- The thin line attached to the 
resonance vertex represents the propagation of a pseudoscalar [g^) or a vector 
meson {gu) and the thick line the propagation of a heavy meson. The outgoing 
line represents a vector meson. 

It is convenient to group together diagrams that are gauge-invariant separately. 
We discuss the various contributions proportional to the coupling constant g^. 
We consider first the terms of class 2) in Fig. |H They are proportional to egp/f 
and involve two types of tensors. We form two gauge-invariant combinations 
with the transverse tensors A°'^''^[p,p) and B°'^''^{p,p), receiving respectively 
contributions from diagrams of classes 1) -|- 2) and 2) -|- 3) of Fig. HI The 
tensors yl°'^''^(p, p) and B°'^'^{p,p) describe the processes where the photon 
couples to the charge of the light or of the heavy pseudoscalar meson, 
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X {Sail + g) (g + 2 ir r (p + 0^ + g^'" (p + 0^} , 
K^fep) = j-^. Sail) S,{p + 1) 

X {StiP + l){p + p + 2l)^' {p + Z)/^ + (^^'^ r] , (88) 

with 

^a(p) = . (89) 

We note that for technical simphcity the tensors ( l88i) are constructed with an 
effective molecule vertex of the generic form D $ ^ as implied by the on-shell 
reduction technique applied in the coupled-channel computation of Section 
2.1. This is the reason why diagrams of class 4) do not contribute to the 
expressions (188|) . The latter are associated with a process where the photon is 
emitted from a hadronic 3-point vertex involving the resonance field. 

The contributions implied by the interactions involving the electromagnetic 
field strength tensor F^^ are proportional to the anomalous coupling strengths 
e^, eg, ec of fl62l [551) or the parameter ep of fl60l) . The terms proportional to 
ep/{f my) probe the class 2) of Fig. Hlonly. Their effect is encoded into the 
transverse tensor A°'^'^{p,p), 



ATip,p) = / Sail) S,ip + /) {(g ■ /) g'^- - P g"} (p + /)^(90) 

Additional tensors C°'^''^{p,p) and C'^^'^{j),p) describe the processes belonging 
to the diagrams of class 1) where the Goldstone boson, which emits the photon, 
is converted into a light vector meson. They are proportional to the parameter 
combinations CAgr/P and CAgE/P respectively. We introduce 



C^fiP.P) = / tSi Sail) S.ip + I) 



<7Tp 

dH 



x2q,l, e^^% r (P + I)' e^./ S^^il + g) , 



C:,'f{p,p) = -^ I 7^,Sa{l)S^{p + l) 



(2vr) 



x2 e'^^^ip + l)p e"%, S^^il + q) , (91) 



where 
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(ml - p^) 



p2 _ ^2 -|- ^ g 

+gf,a Pu Pl3 - 5'm/3 Pu Pa-ifJ'^ l^) 
Sil'^" (P) = Pa Sr^" (p) , S:^''' (p) = Sr^'' (p) pp . 



(92) 



Anomalous processes analogous to those described by the tensors (!9T!) are 
driven by the parameters cq and ec- In this case the photon is emitted by 
a pseudoscalar D-meson, which is converted into a vector D-meson. The con- 
tribution is included in class 3) of Fig. |H We define the corresponding loop 
tensor 



DTiP.P) j^A Sail) Scip + I) 



(93) 



: 2 (/ + p)^ e^^%{/, 6"^% SriP + 1)-P'' I' e,,/ S^'^P + /) } • 



We emphasize that each of the tensor integrals introduced in (!88| [90l [9T| [931) is 
gauge-invariant separately, i.e. vanishes identically if contracted with g^. This 
was checked by explicit calculations. 

Collecting all contributions from the KD and rj channels, we arrive at the 
following decay amplitude. 



"-fhTj^[Al%^{p,p) + Bl%^{p,p)} 



^ hT^ Aft'{p,p) - 2 5;^f (;.,p) 



48 Pmy 

eA 



+ 



+ 



^KD {dT Ck^K*d{P^ P) + 9E CfK^D^P^ P) } 

h^o! [dT C^4.D, {p, p) + Qe C'^4,d^ (p, p) } 



9p 



6 

where exact isospin symmetry is assumed with the coupling constant 



(94) 



h 



(0+) 
KD 



(0+) 



- M2 - m^^ 



Ml 



2/ 
2^3/ 



rat 



9r 



(95) 
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The values of the couphng constants g^^^ and g\^j^^ are given in Table [2j 
Note that the tensors introduced in flHHj [90l [9T| [931) are not antisymmetric 
in the indices a ^ as of notational convenience. While deriving the decay 
parameter (io+_»^i+ according to flS^ . the wave- function projects onto the re- 
levant component, the application of the projection formula (1571) requiring an 
explicit antisymmetrization. 

The contributions of the (j) D* and K*D channels to the decay parameter as 
well as explicit results for the Passareno-Veltman reduction of the tensor inte- 
grals of (l94l) are detailed in Appendix D. Such contributions are proportional 
to Qh- 



5 Radiative decay of axial-vector molecules 



We derive expressions for the radiative decay amplitudes of the -D*^(2460)- 
meson corresponding to the processes 1^ 7 0~, 1"^ ^ 7O+ and 1+^7 1~. 
We follow the same procedure as in the previous section for the D*o(2317) 
radiative decay, adjusting the tensor forms to the hadronic initial and final 
states. 



5.1 (2460) ^7^,(1968) 



The partial decay width Fi+^^q- is expressed in terms of the transition am- 
plitudes M^'l'^ Q_ determined by one number di+_,^Q- as 



= rfi+_,o- et (g. A,) [g^^ (p . q) - p>^ q-} e"^{p, A,) , (96) 

implying 



J- 1+^70 



12n \ 2Mi+ J ' 

with = {p — qy = Mq_ and = M'^+. Like in (187|) we exploit the fact that 
the decay amplitude is antisymmetric in a ^ p. 
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In this case again we have to consider the 5 classes of contributions depicted 
in Fig. m Any diagram is proportional to one of the three coupling constants 
Qji, qh or qh- In the first case, the thin and thick lines attached to the resonance 
vertex stand for the propagation of Goldstone bosons and vector D-mesons. 
In the second case, the role of thin and thick lines is interchanged in the sense 
that the thin lines correspond to light vector mesons, whereas the thick lines 
describe heavy pseudoscalar mesons. In the third case both lines represent 
vector mesons. The outgoing line is always a pseudoscalar Dg-meson. 

We discuss the various contributions, starting with the diagrams proportional 
to the coupling constant qr. There are two types of gauge-invariant combina- 
tions proportional to egp/f. We form two corresponding transverse tensors 
^M,"/3(p^p) and They receive contributions from classes 1) + 2) + 

4) and 2) + 3) + 4) of Fig. H] and describe the processes where the photon 
couples to the charge of the light and heavy meson respectively, 

+p-p Sf{p + /) g''^ + Sa{l + q) p-p S!\p + /) (g + 2 

+P^ Sr^P + Sjip + l)p- - l-p Sf^ip + /) 
+Pp st'^'iP + /) + Pp S!^{p + /) g>^^^ . (98) 

The notations are the same as in (ISUl WM- The tensors fl98p are constructed 
with an effective molecule vertex of the generic form (daD""^) ^ [d^ R/s-r) as 
implied by the projector technique of Section 2.2. The latter is equivalent to 
the term proportional to gn in fl7i|) (see the discussion following fl75l) ). 

The contributions driven by interactions involving the electromagnetic field 
strength tensor F^j, are encoded into the transverse tensor A'^''^^'{p,p), 



Kir'iP,P) = +^ I ^5a(0|(^9)^?^.-/^g.}p/35^^(p + 0p"•(99) 

Additional tensors C'^''^^{p,p) and C^'°'^{p^p) describe the processes where 
the Goldstone boson, which emits the photon, is converted into a light vector 
meson. They are proportional to the parameter combinations CAgrl P and 
^AgEl P respectively. We introduce 
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X e^^^^e^^^^qn^. (100) 

The terms proportional to ec,eQ are described by the tensor D^'"^(p, p). In 
the convention of Fig. H] they correspond to contributions of class 3), which 
probe the anomalous magnetic moments of the vector D-mesons. We introduce 



-Pa SriP + /) qr Sfip + /) . (101) 

Each of the tensor integrals introduced in fl98H10ip is gauge-invariant sepa- 
rately. 

Collecting all terms for the KD* and f] D* channels, we arrive at the decay 
amplitude 



hiil ATn'4P,p) -^"-f h^^} (p,p) 
/ "^v V -J / 



eA 



'Aof2 ^KD* {dT C^'k'D* iP. P) + 9E C^Kf*D' iP. P) } 

__gp_ r 2 3 eg + eg - 3 e ^t,,a(3( ^ 



6ec - en - 3e 



hitDTliP.v)}. (102) 



where 



h^^+} = 72 ^ ^ ^ ^ ^ . (103) 
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The values of the couphng constants g\^j] and g\^jj^ are given in Table [2j We 
reemphasize that the application of the projection formula fl971) requires an 
antisymmetrization of the tensors (198111011) . 

We turn to the contributions implied by the channels involving the K* or 
(p meson. Such terms are proportional to the coupling constants qh or gii 
introduced in (IHOj) . 



While the terms proportional to gn are deferred to Appendix E, we detail 
those proportional to gn here. The former encode the physics of the K*D* 
and (j) D* channels. Appendix E provides in addition explicit results for the 
Passareno- Veltman reduction of the tensor integrals (198111011 11041) . The terms 
proportional to gn are evaluated like those proportional to g^. Formally the 
role of light and heavy intermediate lines in Fig. |4] is interchanged but the 
result involves the same tensor integrals as (198111011) . Additional tensors are 
required to describe the effects of the terms proportional to 7^ or 7^ 
introduced in (1661) . The latter give rise to contributions of class 2) in Fig. |H 
They are analogous to the contributions proportional to ep in (I1Q2I) . We form 
the gauge-invariant tensors p) and B^'"^{p,p), 



Kf{P. P) = -^'J^A ^"(0 S^'^ip + I) p- ip-l). (104) 

We collect the contributions of the K*D and (pD^ channels to the decay am- 
plitude 



egy gn 
P 



2 B^d'^^{j>,p) 



+ 

+ 



[ec + eq/Q) gn 



[ec - eQ/3)gH f ^f,,ai3 „n , ~ ^m,c./3 

g n2j^2 \9eCd^u,^{p,p} + gTCu^D.^{p,p}j 



(105) 



in terms of the tensor integrals of (198111001 1104p . There is no term involving the 
tensor (llOip since we neglect the effect of the anomalous magnetic moment of 
the light vector mesons. 
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5.2 D:i(2460)-^7^:o(2317) 



The partial decay width ri+^^o+ is expressed in terms of the transition am- 
phtudes M{'4.'^^Q+ determined by the number di+^^Q+ as 



di+^,o+ej.(g,A,) qre^^^v" 4^e"^(p,A,), (106) 

VP 



implying 



r 



1+^70+ 



127r V 2 Ml 



+ 



cil+^70+ = (107) 
with = {p — qY = Mg+ and = Ml+. 

There are 4 classes of contributions. They are depicted in Fig. [71 Like in 
Fig. m the solid lines stand for the propagation of pseudoscalar or vector 
mesons. The thick lines are used for the heavy mesons, the thin ones for the 
light mesons. The solid and dashed double lines represent the molecule of the 
initial and final state respectively. Any diagram is proportional to one of the 
four products of coupling constants gngR, gngn, gugn. or gngH- We discuss 
the four possibilities case by case. 

In the first two cases (terms proportional to gR cjr or g^ gn), the thin and thick 
lines attached to the final molecule vertex represent the propagation of pseu- 
doscalar light and heavy mesons. There is only one generic tensor y4^^'"^^(p, p) 
describing these processes. They are of class 1) or 3) in the convention of Fig. 
[71 The contributions probe anomalous electromagnetic vertices proportional 






Fig. 7. Diagrams contributing to the process 1"*" — > 76"^. 
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to ca or 6(7 and eq. We define 



AZ'M p) = -2t qr e^^ .J^, S^) 3,(1 + p) (l + PY Sf{l +p)p'^. 

(108) 

For the terms proportional to gnQRi there is a class 1) contribution only, with 
thick lines describing the propagation of vector mesons. The thin line changes 
from a pseudoscalar to a vector line at the photon vertex. The contributions 
are proportional to the anomalous coupling constant ca- We introduce the 
corresponding gauge-invariant tensor 

DttciP^p) = -2tq.e>^\^g,^l -^J p"" S^) ST {I + q) Sf{p + I) . 

(109) 

For the contributions proportional to gndH, there are two types of gauge- 
invariant combinations. We form two transverse tensors B'^°'^{p, p) and C^'°^(p, p) 
receiving contributions from classes 2) -|- 3) -|- 4) and 1) -|- 2) -|- 4) of Fig. [7] 
respectively. The tensors B^'°''^{p,p) and C^'"^{p,p) describe the processes 
where the photon couples to the charge of the heavy and light mesons. We 
introduce the gauge-invariant tensors 

+sr{p+i)sr'\p+i)]Y 
cz'ip^p) = +^e^^/^7p-p / sr\i) ^p-sni+p) 

+ gf^-Sril+p)+p''Sr^il+p) 

+ g-nnP"- [sf{p + /) srip + + sriP + /) s^^ip + /)}| . 

(110) 

The tensor integrals introduced in (11 08111 lOl) are gauge-invariant separately. 
Collecting all terms we arrive at the decay amplitude 

_ ec + Cq/Q (0+) (1+) , . 

^^^l+^-yO+ — 2M^ "'KD "'KD* ^+,KDD*\P^ P) 
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+ 



ec - eql?, (o+) , (i+) 



+ 



2 My "^^^ '^^^ ^ + ,VDs^s 



{^H ^Sd ^+'!dkk* iP, P) + h^KD D^j<'Kd{p. p) } 



24 P my 



/3 



4/2 

9h (ec + cq/G 



'H9H f ec + eQ/<3 p ec-eQi6 p 



where the couphng constants h^o and /i^^)^ are specified in (195 ^ 11031) in terms 
of the values given in Table [21 Appendix F provides the results of a Passareno- 
Veltman reduction of the tensor integrals fllOStillOj) . 



5.3 L>*i(2460) 7L>*(2112) 



This radiative decay mode is described by a rank-five transition tensor M"/^^"j'!_ 
The decay amplitude has a slightly more complicated structure than those dis- 
played in ( !85l 11061) . It is characterized by two scalar decay parameters as the 
decay may go via an s-wave or a d-wave transition. Separating these two tran- 
sitions is an original and interesting feature of our work leading to angular 
distributions. We write 



d' 



(1) 



(g -p) 



i<l-p), 



Vp^ 



;il2) 



To verify that there are indeed only two independent decay parameters re- 
quires using the identity 



(113) 



We obtain for the decay width 



12 TT 



"l+^7l- 

Mi+ 



+ 



"l+->7l- 

Mi- 



^M2+ -M2_ 
2Mi+ 



(114) 
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with = Ml- and = . As indicated by flllSp , it is quite cumbersome 
to extract the two decay parameters from a given amphtude. Fortunately, this 
task can be streamhned considerably by the projection identities, 



7(1) 



-il6Mi + 



Mi-{Ml+-Ml_f 
(2) _ +i 16 Ml- 

Mi+ (M2+ -M2_)3 



P H ^arafj. Pp -^^^i+^^i- Ha Pf3 , 



P q ^ar/ia Qa Pfi ^^h + ^-y 1' PP > 



;ii5) 



resulting from the antisymmetry and transversality of the decay amplitude. 

We consider the 5 classes of contributions depicted in Fig. HI Any diagram 
is proportional to one of the three coupling constants gn^gn or gu. In the 
first case the thin and thick lines attached to the resonance vertex stand 
for the propagation of Goldstone bosons and vector D-mesons. In the second 
case the thin lines correspond to light vector mesons whereas the thick lines 
describe heavy pseudoscalar mesons. In the third case both lines stand for 
vector mesons. The outgoing line is always a vector D*-meson. 

We discuss the various contributions starting with diagrams proportional to 
the coupling constant g^. There are two types of gauge-invariant combinations 
proportional to egp/f. We form the two corresponding transverse tensors 

, receiving contributions from classes 1) + 2) + 
4) + 5) and 2) + 3) + 4) + 5) of Fig. |4] respectively. The tensors p) 
and B°f'^''^^{p^p) describe the processes where the photon couples to the 
charge of the light and heavy mesons. 



Afr'ip.p) 



+ r Sl^ip + /) (7^° + (/ + qr Sail + q) Sl'ip + l){q + 2 lY p' 



g^''Sf{p + l)p'^ 



+ ip g^^ s^'^ip + /) + r ip sr^iP + g'"' 
+ p^{l + q), Sail + q) s^'^ip + /) (g + 2 IYp^ 
dH 



g%p''sr^{p + l)p'^ 



Sa{l){e^^ 



+ Sl^ip + I) (7^" + F Sl'^'^ip + I) p 



+ LSr^P + l)p'' 



+ r g,^ {si'^'^ip + /) sfip + /) + srip + 1) sr^P + 0} p° 

p, 



+ e^Kr Sf r Sl'ip + l)p-T \ + I, g^^ S^'ip + /) P° 
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+ IpP^ Sr^P + I) 9''^ + IpP^ S^'^P + I) 

+ IpP"" 9-nn [Sr%P + /) Sr^ip + l) + Sr^^P + /) Sf{p + l)}p- 

(116) 

In the tensors A^'^^j^'"^ and B'^^^^'"^ there are contributions where the photon 
couples to the final vector particle. In the convention of Fig. Ill these are 
diagrams of class 5). They give rise to terms proportional to Sp^'^^^^p), with 
rrip = p^. Such contributions are not at odds with parity conservation since 
the tensor field carries spin one quanta with both parities. Analogous terms 
where the photon couples to the initial vector meson do not arise due to parity 
conservation. Since the resonance field couples always with [drR^"'), only the 
positive parity component is accessible. 

The contributions induced by interactions involving the electromagnetic field 
strength tensor F^^, are proportional to the anomalous coupling strengths e^, 
eg, 6(7, eg and ec of fl^ ISSl IF?]) or the parameter cp of (pPj) . The terms 
proportional to ep/{fmy) probe the class 2) of Fig. |4] only. Their effect is 
encoded into the transverse tensor A'^^''^''^^'{p,p). The tensor B^j^"'^'"'^ {p,p) is 
associated with contributions probing the anomalous magnetic moment of the 
D mesons, which are proportional to ec and eg. 



T e.. (l" q,-{l- q) 9"^ ^ S^^'^^p + I)] p\ 



(117) 

B:t''^P.v) = +^ J [l''q.{sr{P + l)Sf{p + l) 

- S7{p + /) Sf{p + /)} p"] - ejP [l, q^p^ {Sr\p + /) Sf{p + /) 



-sr'^{p+i)sf{p+i)]p'' 



:ii8) 



The additional tensors C+^'^'"^(p, p) and C°'^'^''^^{p,p) describe the processes 
where the Goldstone boson, which emits the photon, is converted into a light 
vector meson as included in class 1). They are proportional to the parameter 
combinations eAgv/P and eAQr/P respectively. We define 



cgai:f'iP, = +2 ^ / ^-^(0 .r \sr{i + q) s!\i + p) f 
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± r 9-nn St {I + q) Sf'/'il + p) j , 

(119) 

xg,,sf'''il + q)Sfil+p)p-. (120) 

Anomalous processes analogous to those described by the tensors flll9p are 
driven by eg and ec- In this case the photon is emitted by a pseudoscalar D- 
meson, which is converted into a vector D-meson. The contribution is included 
in class 3) of Fig. |H We introduce the corresponding loop tensor 

X r / e'^- {l+pr Sf{l +p)p^. (121) 

The tensor integrals introduced in (11 16111 171 II 191 fT2Ti) are again gauge- invariant 
separately. 

Summing all terms, we get the decay amplitude 

+47^''™'-^*- (P'P' + ^^"K'^Ssf (ftp) 

- T^T^ h^Kt [gv cftK^n^ ip, v) + gr cfirJt iP, P) } 

+ '-^^h^^llDf--fip,p)}. (122) 

We turn to the contributions induced by the coupling constants gn and gn 
introduced in (IHOl) . The terms proportional to g^ are deferred to Appendix 
G. We discuss those proportional to gn- The former are expected to be more 
relevant than the latter due to phase space. Appendix G provides in addition 
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explicit results for the Passareno-Veltman reduction of the tensor integrals 

(IIMmlimD. 



The evaluation of terms proportional to ^ji^ is similar to those proportional 
to Qff. Formally the role of light and heavy intermediate lines in Fig. H] is 
interchanged so that the tensor integrals formed in flll6H121]) will occur again. 
Additional tensors are required to describe the implications of ct 7^ as 
introduced in fl66l) . The latter give rise to contributions of class 2) in Fig. |H 
They are analogous to the contributions proportional to ep in (11111) . We form 
the gauge-invariant tensor B^''^^{p,p) 



B 



ab 



dH 



p 



(123) 



The contributions proportional to gn are 



_ e{gE-gT)~gH ^) + B^_%r^^ip,p)} 

-i^Bf--\pM 

-'~^''{igE + -gr) Al%^fip,p) + {gE - -gr) fep)} 



A f 2 I ^'^ ' '-'+,DD*K* [P^ P) + ^gv '^DD*K* IP' P) 



+ {gv-gT)C^.%'b^.Up.p)] 
+Cgv-gT)c^-%rDUp,p)} 

+ ^^f^ ^^^r/fe P) + "^ff^ Dt-r^P. P) . (124) 



6 Numerical results 



We confront the results of the previous sections with the experimental data 
on the radiative and strong decays of the scalar and axial-vector Dg-mesons 
reviewed in the introduction. 
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Our prediction of 140 keV for the isospin-violating strong width D*q(2317)^ — > 
Ds(1968)^ vr" is compatible with the empirical bound F < 3.8 MeV [23] but 
this comparison does not provide any significant constraint on the underly- 
ing coupled-channel dynamics. The present upper limit on the ratio of the 
radiative to pionic decay width [21] 



r[g; (231^^ - "tt < 0^059 (125) 

r[D*o(2317)± ^ D,(1968)±7rO^ ^ ' 



implies for the decay constant (io+-+7i- defined in f lS^ the inequality, 



0+^71- 



< 0.117 GeV^S (126) 



using our predicted value for the Z)*q tt^ decay width (140 keV). 

The total width of the 0*^(2460)=^ meson is less than 3.5 MeV [23]. Our pre- 
diction of 140 keV for the isospin-violating strong decay width D*i[2460) 
D*(2112) vr" is compatible with that upper bound but this is again not very 
significant. The constraints on the radiative decays of the D*^(2460)^ to the 
D,(1968)±, the 0^(2112)^ and the 0^0(2317)^ given in (ISD, 1^ and ([5]) imply 
the following relations for the decay constants di+^^Q-, di+_^Q+ and d\+_^^^- 
introduced in ([97]), ffTOTD and ffTTill . 



(iiH 



>7 0- 



0.138 GeV-1 



d 



1+-+70+ 



< 0.665 GeV"^ 



"l+-^7l- 



+ (1.164)^ 



"1+^71- 



< 0.39 



(127) 



It is interesting to confront the empirical constraints (11261 11271) to the predic- 
tions of heavy-quark symmetry. To leading order, the heavy-quark symmetry 
implies the relations 



.7I' 



»7 0- 



(1) 

1+^71- 



d 



(2) 

1+^71- 



d . 



where Mc — My is a typical mass of a charmed meson. The result ( 1128P can 
be derived from [10] by a matching of corresponding decay amplitudes. In the 
heavy-quark mass limit the parameter d scales linearly with the charm quark 



mass. The value given in fll27p for 



0.30 for Mr 



diA 



.70- 



suggests the range 0.25 < \d\ < 



2000 MeV. This value is barely compatible with the 



bound (I126P and requires that 



di 



0+-+71- 



be close to that upper bound. 



We discuss the electromagnetic decay parameters obtained using the expres- 
sions derived in Sections 4 and 5. The contributions from the KD and r]Ds 
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channels are shown in Table [3] for radiative transitions between scalar and 
vector states and in Table H] for the radiative decay of the 1^ to the 1~ state. 
The first column displays contributions involving the 4-point vertices (1541) 
and (1601) . The second column shows the effect of anomalous processes and the 
third column contributions induced by the light vector mesons as intermediate 
states. 



We consider the first column of these tables and try to reproduce the con- 





10 X do+^^i-[GeV-^] 


KD 


+2.371 5P + 3.036 ep 


-2.252 (ec+eQ/6)5P 


-{6.095 gr + 0.950 qe) CA 


vDs 


+6.822 5p 


-1.018 (ec - eQ/3) gp 


-iU.OSgr + 3.575 gE)eA 




10 X di+_^o-[GeV-i] 


KD 


+1.641 5p + 2.920 ep 


+2.272 (ec-eQ/6)(7P 


+ (0.500 + 0.501 g£;)eA 


r]Ds 


+3.445 5p 


+1.222 (ec + eQ/3) gp 


+(3.1495T + 0.883 5i?)eA 




10 X di+^^o+[GeV-i] 


KD 


+2.104 (ec + eQ/6) 




+0.891 {ec - eQ/3) 


Table 3 

Contributions to the decay constants dg+^^i- and di+^ 
and (|111|). We use the coupling; constants of Table [2l f = 
and Mv = 2000 MeV. 


^0+ implied by ([Ml), (|102D 
= 90 MeV, mv = 776 MeV 






KD 


+5.U7gp + 7.347ep 


+0.221 (ec + eQ/6) gp 
+3.214 (ec -eQ/6)5P 


+ (8.544 - l.i58 gv)eA 




+11.92 ffp 


-0.227 (ec- eQ/3) gp 
+1.756 (ec + eQ/3)5P 


+{26.85 gT-5.9Ugv)eA 




10 X df}_^^^_ 


KD 


+2. 799 5P + 6.305 ep 


+2.153 (ec + eQ/6)5P 
+0.0650 (ec - eQ/6) gp 


-{8.802 gr - 1.335 gv)eA 


vDs 


+8.852 5p 


+0.489 (ec - eQ/3) gp 
+0.080 (ec + eQ/3)5P 


-{21.62 gT- 3.469 gv)eA 



Table 4 

Contributions to the decay parameters i- and -|^_ implied by (|122p . We 

use the same coupling constants as in Tablefe] 
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straints fll26p and (11270 using these terms only and the empirical value gp — 
0.57. This exercise leads immediately to the conclusion that gp and cp must 
have opposite signs. Adding constructively the contributions from the KD 
and 7] Dg channels would produce decay parameters that are much too large. 
We recall that ep parameterizes gauge- invariant 4-point vertices describing 
the process where a D-meson emits a photon and a charged Goldstone boson 
simultaneously. Such interactions must be taken into account in effective field 
theories. The parameter ep can be varied to achieve consistency with (11261) 
and (I127p . To fit the 1+ ^ 7 0~ decay, the range of values for ep is limited to 
the intervals —0.57 < ep < —0.49 and —1.50 < ep < —1.42. The constraints 
on the process 1+^71" suggest —1.54 < ep < —0.84 for gp = gp = 0.57. 
From the 0"*" decay we deduce —2.08 < ep < —1.37. With ep = ep ^ —1.50 
we arrive at an acceptable scenario, given the assumption of Section 3 that 
the chiral power assigned to the ep vertex is promoted from order Q^, to or- 
der Or,. The corresponding decay parameters are collected in the first column 
of Table [51 We observe that this scenario, though compatible with the con- 
straints (I126p and (11270 . is characterized by decay parameters in significant 
disagreement with the heavy-quark symmetry relation (11280 . 

A value of ep of the order of —1.5 has many consequences. As mentioned in the 
discussion following (16T|) . this coupling is unnaturally large and points to the 
effectiveness of the theory. It also induces very important cancellations in the 
KD channel and between the KD and the rjDs channels. The 1^*0(2317) — > 
7 (21 12) transition is actually dominated by the contribution of the rjDg 
channel. The corresponding radiative width is 2.85 keV and the calculated 
value for the ratio ([2]) is 0.02. The comparison with other coupled-channel 
calculations is not easy as the mechanism driving the D*q{2317) 7 D*(2112) 
transition is not simple. Our radiative width is larger than the value of 0.49 keV 
obtained in [TB] but we note that similar destructive effects between channels 
are observed in both approaches. Our result is quantitatively comparable to 
the radiative width of ~1 keV found in [12] in a molecular picture, despite 
the fact that the r]Ds channel is not taken into account in [19] and important 
in our work as well as in [18]. This particular example illustrates the need 
for a consistent description of the spectroscopic properties of the D^-mesons 
guided by general principles such as strong interaction symmetries. For the 
decays ^^1(2460) 7^,(1968), Dl^{2AQ0) ^ 7L>,o(2317) and D*^(2460) 
7 (21 12), we obtain partial widths of 50 keV, and 18 keV respectively 
(assuming always a strong width of 140 MeV). We stress that performing a 
formal expansion of the full expressions for the decay parameter in the inverse 
charm-quark mass leads to results that are compatible with the expectation 
from the heavy-quark symmetry. With phenomenological charm quark masses, 
the decay parameters obtained for the KD channel distort the pattern implied 
by the relation (11280 obtained in the limit of infinitely heavy charm quarks. 
Such breaking pattern is even larger in the rjDs channel. We note also from 
the numbers provided in Tables [3] and H] that the contribution proportional to 
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ij LL) 111 J J- * J 


"0+^71- l^-'OV J 


_i_n nRQ _i_n (t?^ n nfr^ -i_n n?'^ 
-|-u.UDy -|-u.uoo — u.uuo -rU.u/o 


"1+^70- I'-^cv J 






+0.055 +0.055 +0.055 


"1+^71- 


-0.129 -0.097 -0.105 -0.090 


"l+-»7l- 


-0.282 -0.251 -0.251 -0.251 



Table 5 

Decay constants that are implied by ep = ep = —1.50 and gn = 9h = dH = 0- 
The values used for e^, ec = ec and eg = eg are discussed in the text. The four 
scenarios are characterized by I) = and eg = ec = 0, II) ca = 0, III) ca > 
and IV) CA < 0. We use / = 90 MeV, my = 776 MeV, My = 2000 MeV. 

gp is larger in the rjDg channel than in the KD channel. This is because there 
is only one contribution in the rjDs channel while there are two contributions 
of opposite sign in the KD channel corresponding to the graphs where the 
photon couples to the or to the D~^. 

We display the effect of the anomalous contributions in the second column of 
Tables [3] and HI Without light vector mesons as explicit degrees of freedom, 
these terms determine entirely the l"*" — > 7 0+ process. According to Appendix 
A, we should take ec — 0.13 and adjust the value of eg to the channel under 
consideration (eg ^ 0.91 for the KD channel and eg ~ 0.52 for the rjDg 
channel). Given the coupling constants of Tabled we predict di+^^Q+ = 0.055 
GeV~^, a result compatible with the bound of f ll27p . 

We discuss now the effects induced by the presence of light vector mesons 
when they do not couple directly to the scalar and axial-vector molecules 
idH = gH = 9H = 0) but influence their radiative decays as intermediate 
states through the anomalous vertex introduced in (162|) . This contribution 
is proportional to ca and given in the third column of Tables [3] and HI The 
parameter ca exhibits significant flavour SU(3) breaking as indicated in fl64p . 
To take it into account we use 

ej^=) = e^*^^^)~ ±0.053, 

ef = 2er°^^^°) -er"^^"""^ ^ ±0-148, (129) 

in the rjDg and KD channels respectively. The phase of the parameter ca is 
not determined by experiment. We consider both signs and provide the nu- 
merical values of the electromagnetic decay parameters in the third and fourth 
columns of Table [5] for e^ >0 and ca <0. The phases and the size of the pa- 
rameters qt, Qei Qt and gy were estimated in Section 3 from the assumption 
of universally coupled light vector mesons, together with the ansatz of a com- 
bined heavy-quark and flavour SU(4) symmetry. We used the values gv — 0.71, 
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dE = Qp and Qt = Qt = 0.5 mv gp/ My with gp = 0.57. We caution that these 
numbers were derived with significant approximations. As can be seen in Table 
the sign of matters mostly for the radiative decay width of the 1^*0(2317): 
this width is extremely small for ca >0 (5.4 eV) and much larger for <0 
(3.2 keV). An accurate measurement of the D*]^(2460) 7D*(2112) decay 
width would therefore be most helpful in determining whether the Lagrangian 
(1621) leads to constructive or destructive interference effects in that quantity. 
We note that the and channels contribute to the decay ampli- 

tudes with opposite signs, the contribution of the K'^ channel being twice 
as large as the contribution of the channel. 

We enlarge the discussion of the role played by light vector mesons by assuming 
that the K*D* and D* channels can couple directly to the scalar and axial- 
vector molecules and the K*D and Dg channels to the axial- vector molecules 
through the interactions defined in (IHOl) . In the heavy-quark mass limit, the 
coupling constants are equal, i.e. gn = gu = gu . 

We consider first the transition 1"*" 7 0"'" between molecules. The corre- 
sponding decay parameter reads 



10 X (ii+_^o+[GeV ^] =gHgH 



1.187- 0.377 



-gH gn 



5.363 (ec + cq/G) + 5.607 (ec - cq/S) 



0A29 gH + 9.2AA cjH 



JvDs) 



1.704 (7// + 13.62 , (130) 



where we separate the contributions from channels involving the D and Dg 
mesons. In the first term, the numbers 1.187 and —0.377 represent the contri- 
butions from the (pD* and the K*D* channels respectively. Similarly the con- 
tributions proportional to ec + gq/G and — eg/S reflect the K*D — * K*D* 
and (pDg —>■ (p D* transitions. With the values of e^, eg and ec discussed 
above, we may use (11301) and (11271) to constrain the parameters gH,gH and 
gn- In the heavy-quark limit {gn = gn = gn), we derive the conditions 



- 10.5 <gH < 2.18 , -2.18 < gn < 10.5 , (131) 

for the positive and negative values of respectively. Clearly the result (11311) 
does not provide a strong constraint on the size of the parameter gn- 

We turn to the radiative transitions of the scalar and axial-vector molecules to 
0~ and 1~ ground states. We display in Tables [H] and [7] the contributions of the 
channels involving light vector mesons to these decays. The terms shown in 
Table [6] correspond to the expressions (I165P and (I105p (obtained in Appendix 
D and in Section 5) for the 0"*" — > 71^ transition and (I169P derived in Appendix 
D for the 1"*" — > 70~ transition. The terms quoted in Table [7] correspond to 
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in X fln4- 1 fOpV~^l 


K*D* 
^ 9H ■ 


+10.66 gr + 2.976 +1.758 er + 2.343 - 2.765 ey 
-27. IS EA gp +(ec + eg/e) (3.408 + 0.611 ge) 

+ (ec7 - eg/bj (2.283 5t + 0.653 g'v'j 


^ n* 

~ 9H ■ 


-60.24 e^^P +{ec - eq/^) (1.383 + 0.322 gg) 

+fg^ + (I 278 + 592 oi/l 




10 X (ii+_^o-[GeV-i] 


A JJ 
~ 9H ■ 


+(ec7 + eQ/6) (6.777 5r + 1.337^^^) 


4>Ds 
~ QH ■ 


-2.902 5y +6.462 eg 

+ {ec — eQ/6) (3.548 fi(T + 0.84/ gE) 


K*D* 
~ QH ■ 


+9.798 + 1-603 gE +3.338 er 

-9.070 e^SP +1.568 (ec7 + eQ/6) ffy 

— (ec — eg/bj (2.(^Uo(7t + U.o4(^ gE) 


<PDI 
~ QH ■ 


+16.92 gr + 3.807 gE 

-15.44 e^SP +1.072 (ec -eQ/3)ffy 

-(ec + eQ/3) (1.478 5T + 0.359 gij) 



Table 6 

Contributions to the decay parameters (Iq+^^i- and di+_,^Q- implied by the light 
vector meson couphngs as given by (jl65l 11051 1169p . The various terms have to be 
multiplied by either gHidH or gn as indicated in the first column. 

the expressions fll24p and fll79p (obtained in Section 5 and in Appendix G) for 
the 1"*" — > 7!" transition. They involve the additional parameters e^, e^, ey, ey 
and eE,eE introduced in (166|) . Clearly, it is not possible to determine these 
numbers by the constraints (11260 and (I127p only. 

To achieve a qualitative understanding of these terms, we assume that all 
parameters are correlated as dictated by the heavy-quark mass limit (i.e. ep = 
GPi^v = ey,...). In this limit, there remain five unknown parameters gn-, ep, 
ey, ct and ce- One may expect to learn little from the constraints (11260 and 
(11270 only. However, this is not quite so. For specific values of qh and ep, 
one can always adjust the values for ey, and ce so as to reproduce a given 
value of the decay parameter (io+-»7i- and the partial decay widths of the 
process l"*" — > 70~ and 1+ — > 7 1~ derived in (11260 and (11270 . We assume as 
before gv — 0.71, qe = gp and gr = gr = 0.5mv gp/My with gp = 0.57. 
The requirement that the parameters ey,eT and eg be real, as implied by 
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1+— +7 1 


A U 
~ QH ■ 


1 1 no/i ;^ 1 n 1 c^o ^ i n ^i c 
+ i.Uz4 (77^ + U.iOZ (7£; +U.Dioer 

+2.049 CA 9P -{ec + eg/e) (2.930 (7t + 0.483 ^y) 




+2.67A gr + 0.597 gE 

+5.7i4eA gp -(ec - eg/dj (i.bOb^T + 0.147 gyj 


K*D* 


+1.065 5r + 0.569 gv +0.528 er - 0.979 ey + 1.069 e^; 
-5.914 CA gp -{ec + eQ/6) (0.209 + 0.037 ffjj) 

-\-{ec — eQ/b) [{)A6\)gT + U.UU^^yj 


^ n* 
~ 9H ■ 


J-9 ^Q7 /^rn -L n 7Af\ n-r r J- 1 "^98 

-t-Z.oy / y'jp -\- \j. 1 yy -\- l.o^o o ^ 

-10.99 eA gp -{ec - eg/S) (0.047 + 0.011 gs) 

-i-ien + enl?,'] fO 260 ot + 058 a^A 






A U 


1 n fiCQ p; 1 n nni i n c:oi p^ 
H-U.Doo 5^7^ + U.Uyi Qe +U.0Z1 Ct 

-6.998 eA9P +(ec + eQ/6) (1.567 5r - 0.202 gy) 


~ 9H ■ 


+2. 276 + 0.5045s 

-13.99eA gp +{ec - eQ/6) (0.793 - 0.190 gyj 


K*D* 
~ 9H ■ 


+1.512 gr + 0.446 gy +0.669 er - 0.907 ey + 0.881 e-E 
+1.435 eA 9P +(ec + eg/e) (0.247 + 0.044 ge) 

— \.ec — eQ/bj — \J.y}6i gv) 




+2.8205T + 0.5875y +1.041 

+2.858 eA ~gp +(ec - eQ/3) (0.1575T + OMT gE) 

-{ec + eQ/3) (0.443 (7T - 0.014gy) 



Table 7 

Contributions to the decay parameters d^i+^^i- and ^[+^71- implied by the light 
vector meson couplings as given by p24l [T75|) . The various terms have to be mul- 
tiplied by either gn or gu as indicated in the first column. 

charge conjugation symmetry, defines stringent conditions on the ranges of 
values allowed for gu and ep. Note that the latter are determined by quadratic 
equations, the solution of which involves a square root. The constants ey, ct 
and ce be real numbers only if the argument of that square root is positive. 
For a given value of ep, we derive the condition that the parameter qh has to 
be confined in a small interval 

crit,— , , crit,+ /looN 

Qh <9h < ■ (132) 
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ij iij 


"0+^71- ['-'CV J 


_i_n nf^7 _i_n 1 n/i 


ai+^^o-[Kje\ J 






-0.043 +0.091 


1+^7 1 


+0.303 +0.278 


1+^7 1 


+0.202 +0.071 



Table 8 

Decay constants obtained using ep = ey = ex = = and ep = ey = ex = = 
0. In set I) we use ba > with gn = 9h = Qh = —0.46. Set II) assumes e^i < 
and QH = gH = 9H = -0.25. We assume / = 90 MeV, my = 776 MeV, My = 2000 
MeV. 

We discriminate between the scenarios impHed by positive or negative values 
of the parameter ca-, for which we find respectively 

^crit,± _ _o,39o + 1.332 do+^71- [GeV^^] - 0.242 ep ± 0.131 , 

^crit,± _ _o.285 + 0.751 rfo+-.7i-[GeV"^] - 0.137 ep ± 0.074 . (133) 

We checked the stability of this result against reasonable variations of the 
parameters Qv^Qt and qe- In (11331) we allow the decay parameter c/g+^^i- 
to take any value. The result (I133p is very significant: for any reasonable 
range of ep it requires the coupling constant qh to be quite small, typically 
\gH\ < 0.5. This justifies in retrospect the coupled-channel computation of 
Section 2, which assumed that the light vector mesons are not relevant for the 
formation of the scalar and axial-vector Dg molecules. 

The light vector mesons may nevertheless change the radiative decay para- 
meters significantly because of cancellations between terms. Our main result 
is that the introduction of light vector mesons leads to a very consistent pic- 
ture. We obtain values for all decay parameters that are compatible with the 
empirical constraints using vanishing values for all gauge-invariant counter 
terms ep = ey = ey = e^ = and ep = ey = ex = ce = 0. In such a sce- 
nario there is one free parameter only, Qh = gn = (Jh, which can be dialed to 
recover all empirical constraints. The results for the positive and negative ca 
scenarios are collected in Table [HI It is interesting to observe that we predict 
a negative decay constant for the l"*" — > 7 0~ decay contradicting the naive 
expectation of heavy-quark symmetry. We emphasize that this follows even 
though performing a formal expansion of the full expressions for the decay 
parameter in the inverse charm-quark mass leads to results compatible with 
the expectation from the heavy-quark symmetry. Such an expansion assumes 
for instance <C Mc, which is not realized too well in nature. Our results 
provide a physical justification for promoting the counter terms proportional 
to ep and ep to chiral order Ql rather than the expected power Ql.. Once the 
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light vector mesons are introduced as important physical degrees of freedom, 
the naturalness assumption for the residual size of ep and cp appears justified. 

We emphasize that it does not appear possible to predict precise values for 
the decay parameters. The results of Table [8] should be viewed as possi- 
ble and natural scenarios. We provide therefore the radiative widths asso- 
ciated with the decay constants of Table [HI assuming the strong width of the 
D*o(2317)=^ and D*i(2460)=^ to be 140 keV, as a mere indication of their ex- 
pected range. For the scenarios labeled I and II, we find ro+_»-yi- = 1.94 and 
6.47 keV, Ti+^^q- =44.50 and 45.14 keV, ri+^^o+=0.13 and 0.59 keV and 
ri+_,^i-=21.8 and 12.47 keV respectively. Precise unquenched lattice QCD 
simulations for the hadronic coupling constants of the Goldstone bosons and 
light vector mesons to the D-mesons would be very helpful. 



7 Summary 

Based on the chiral Lagrangian properties of scalar and axial-vector meson 
molecules with open-charm content were studied. Chiral correction terms were 
incorporated systematically in the coupled-channel dynamics, where we relied 
on constraints from large- iVc QCD and the heavy-quark symmetry. We focused 
on the 1^*0(2317) and -D*^(2460) states and computed their isospin-violating 
hadronic and electromagnetic decay widths. In order to establish manifestly 
gauge-invariant results for the electromagnetic decay parameters the spin-one 
particles were represented in terms of antisymmetric tensor fields rather than 
by the more conventional vector fields. The role of explicit light vector mesons 
in the radiative decays was investigated. 

The main findings of this work for the strong and radiative decay widths of 
the D*o(2317) and Z?*]^(2460) states can be summarized as follows. 

The hadronic isospin-violating decay widths of both the D*o(2317) and i5*i(2460) 
states are predicted to be 140 keV. Chiral corrections to order lead to 
a significant enhancement of these widths. We point out the importance of 
treating consistently the 77 vr mixing and the isospin-mixing effects in the KD 
system. The observed upper limits on these strong widths of 3.8 and 3.5 MeV 
respectively do not provide any useful constraint on the theory yet. 

The radiative decay parameters of the D*o(2317) ^7/}* and L)*^(2460) 
7/^5,7-0* and D*^(2460) 7L'*q(2317) were computed in the hadrogenesis 
conjecture. They involve only one-loop diagrams. We find that both the rjDs 
and ?7-D* channels contribute significantly and that large cancellations occur 
between the different channels. The results are compatible with all empiri- 
cal constraints if one gauge-invariant contact term is considered to be more 
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important than expected from a naive naturalness assumption. The decay pa- 
rameters are subject to a cancellation mechanism, which makes a prediction 
of their precise values unreliable at the moment. 

When considering light vector mesons as explicit degrees of freedom coupled 
to the scalar and axial-vector molecules, radiative decay parameters that are 
compatible with the empirical constraints, can be obtained without envok- 
ing subleading contact terms. We understand this effect as an indication of 
the dynamical role of light vector mesons in the electromagnetic transition 
processes. We predict that the DIq{2317) and -0*^(2460) have small coupling 
strength to the K*D* and (p D* channels despite the decisive role played by 
these degrees of freedom in the radiative decay processes. 

Finally we point out that the invariant t] D* invariant mass distribution shows 
a signal of a member of an exotic axial- vector sextet at a mass of 2568 MeV and 
with a width of 18 MeV. While that state decouples from the tiD* spectrum, 
its heavy-quark partner defines a narrow dip at a mass of 2410 MeV and with 
a width of 2 MeV in the tt D mass distribution. 
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Appendix A 



The modulus of the coupling constant gp appearing in the Lagrangian (1371) 
can be determined from the measured strong D*-meson decay widths [21] 

Td*+^d^^+ = (65.0 ± 14.9) keV , 

Td*+^d+ nO = (29.5 ± 6.8) keV . (134) 
Using 

= 39.60 MeV, 
g^m = 38.30 MeV, (135) 



\9p\ Qcm 



12 vr P 

Ifi'Pp Qcm 



24 TT f 
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leads to 



|c/p| = 0.57 ±0.07. (136) 

The modulus of the coupling constants eg and cc introduced in the Lagrangian 
(155!) can be determined from the observed branching ratios for the radiative 
decays D*+ D+ 7, ^ 7 and D*+ ^ 7 [21] . The total widths of 
the and of the D*~^ are not known. They are the sums of the strong and 
radiative widths. The strong widths associated with 7r° emission, 



I 123 

rz,*o^,,o^o = 1^ % ^ 42keV, gem = 43.12MeV, 
24 vr 

r^:+^z)+.o = ^|^^ ^e2 77keV, g,^ = 47.87 MeV , (137) 

are calculated using the value of \gp\ derived above, e is the ttqI] mixing an- 
gle (fTOl) . From the data on the D*~^ D^-f partial width and the known 
branching ratios of the and D*~^ into the ir^D^ and n^D^ channels 
we obtain 



Td*+^d+^ = (1.5±0.8)keV, 

Td*o^do^ = {TD*o^DO'y + 42keV) x (0.381 ± 0.029) 

= (26.0±3)keV, 
Tnt^^Dt, = i^Dt^-,Dt, + 77keV) x (0.942 ± 0.007) 

= (1328 ± 180) keV. (138) 

From the Lagrangian fl5^ . we derive 



Ml^ ieQ-3e cf 
47r 9Mt 



r^.+^^+^ = ' ^^^"^ q'^^ , q,m = 135.78 MeV , 



M2,o47r(2eQ + 3ec 



i2 



r^.o^^o, = 1cm , Qcm = 137.16 MeV 



9My 



^Dt^^Dt, = ^ ^""^g^r^' qL , Qcm = 138.91 MeV . (139) 



From the first two constraints of (11381) . we get 



leg - 3 eel ^ 0.537+°;J28 , |2 eg + 3 ec| ^ ^-^Ol^oUt , 
i.e. eg = 0.91 ±0.10, ec = 0.13 ± 0.05 , (140) 
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where we assume My = 2000 MeV. Using flMOp and taking e ~ 0.01 in the 
third constraint overestimates the D*~^ — >■ Df'j decay rate by an order of 
magnitude (~1.9 keV rather than ~ 0.13 keV). This may be a signal of SU(3) 
breaking effects in eg or suggest the need for an increased value of e. For 
e = 0.01 and ec = 0.13, an effective eg ~ 0.52 reproduces the empirical 
width. 

The modulus of the coupling constant ca appearing in fl62l) is determined from 
the radiative decay of the light vector mesons, K*^ K^li K*^ K^j and 
(f) ?77. From the Lagrangian fl62l) . we obtain 



i/f.±^A'±7 9/1 ^. 1/1/1 ^ /2 ^2 



24 X 144 TT mf /2 \^ m^,± 



24 X 3671 m^yP \ m\,, 

I 19 5 / 2 \ 3 



^ 24 X 54 vr /2 I ml 



fl411 



where we will assume / = 90 MeV and my = 776 MeV. The experimental 
decay widths [21] imply conflicting values for ca- We find 



Tk*o^ko^ = (116 ± 10)keV, 
TK,±_^K±'y = (50 ± 5) keV , 
r,^^^^ = (55.38 ± 1.68) keV, 



\eA\ = 0.119 ±0.006, 
|e^| = 0.090 ±0.004, 
e^l = 0.053 ±0.001. (142) 



Appendix B 



To work out the implications of the heavy-quark symmetry of QCD, it is useful 
to introduce auxiliary fields, P±{x) and P±{x) such that 



D{x) = e"' P+ (x) ± e+' ('""^ P_ (x) , 

D'^'ix) = i e-^(''-^)*^- {v" P^ix) - P^(x) ± — (a^'P^ - S'^P^) } 

^■^+i{vx)M, {t;'^P^(x) - v" P^{x) - (S'^P^ - S'^P^)} (143) 

with a 4- velocity normalized by t;^ = 1. The mass Mc is an averaged value of 
the pseudoscalar and vector charmed meson ground state masses. The fields 
P±{x) and P±{x) are therefore varying slowly in space and time and their 
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derivatives are small compared to M^Va P- In the limit Mc — ^ oo, these deriva- 
tives can be neglected. 



We rewrite the interaction terms introduced in (1371 [60l [65]) in terms of the 
auxiliary fields. Restricting the Lagrangian density to the 'plus' components 
and dropping that index for simplicity, we have 



C = I [P, (9^$) P-P (9^$) p4 + I e^^'^^ P, [d^^) Pp 

- ^ V, P {dj^-^) P + ^ V, P, {d^V^n 
+ 1 e^^"^ {P^V,,P-P V,, Pp} +t P, V^' P, 

+ i^M^ {p^ Q] P - P [(9.$), Q] P^} 

+ ■■■, (144) 

where the ellipses stand for additional terms involving derivatives of the soft 
fields P = P+ and P^ = Pl^ and the 'minus' components. In deriving (11441) 
we used the equation of motion of the auxiliary field P!^{x), 



TtM,v^P^ + d^P^ = 0, (145) 
obtained from the equation of motion of the field D^", 

df'daD"'' - d'^dc.D"^ + Ml Df"" = 0. (146) 

The QCD action depends linearly on the charm-quark mass. The effective La- 
grangian fll44p should therefore be proportional to M^. The coupling constants 
gpy and gpy must then approach a finite value in the limit of infinite charm- 
quark mass. The same argument requires the coupling constants qt and qt to 
scale with l/Mc- The fact that there is no term proportional to qe displayed 
in (11441) signals that the corresponding interaction term is at least subleading 
in the heavy-quark mass expansion. Note that this observation is compatible 
with a possible finite and non-zero value of the asymptotic value oi qe in the 
heavy- quark mass limit. 

In the limit of infinite quark mass, the fields P± and P± may be combined into 
a multiplet involving the 1~ and 0~ fields. Therefore the properties of pseu- 
doscalar and vector states should be closely related. We follow the formalism 
developed in [2B|27ll28ll29f3Uj and introduce the multiplet field H, 
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H = ^oHho= (4,M 7^ + 4 ^ 75) ^ (l + 1^) , 

Pltv, = 0, v' = l, (147) 

in terms of which the interaction Lagrangian can be recast. Let S denote a 
transformation of the heavy-quark spin symmetry group SU(2)^, whose ele- 
ments are characterized by the 4-vector 6*" with 6 ■ v = 0. According to 
under such a transformation the field H transforms as 



^" = ^75[|^,7l, 51 70 = 70 5a, [^,5J_ = 0. (148) 

Under a Lorentz transformation characterized by the antisymmetric tensor 
ijj^y, the spinor part of the fields transforms as 



^ _^ S^, uji^- ^ g-i S^, _^ o;''- ^ g-i 5^. a;"" 



5;.. = ^[7., 7.]- (149) 

Only the combinations where Dirac matrices are to the right of the field H 
or to the left of the field H are invariant under the spin group SUy{2). It is 
straightforward to construct those terms that can be matched to the structures 
given in f ll44p . We introduce 



£ = -^tr (i7(9^$)75 7^iy 
-^tr {H{d^V-n^,H 



^yi^^'tr (if 75 7^(5.$), Q 



H 
H 



h 



^i^F^^tT [HiVf"", Q}H 



(150) 



where we note that the field if is a three-dimensional row in flavour space, 
each of its components consisting of a 4-dimensional Dirac matrix. We have 



tr75 7/^7i^7Q7/3 = -4^e^l.Q/3 , (151) 
in our convention. Matching the expressions (11441 11501) we obtain 
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yp ^^J-c 


yp ^^J-c 


= fp, 


ep ivic 


Pr, M 


f 

J 


f 

J 


f 77?,?, 


f 77?,?, 


riir M 

yv ^^J-c 


_ yv ^^^c 


= fv, 




Ptr M 


f 

J 


f 

J 


f 777?, 

J '"'V 


f Tn?r 
J '"'V 




grM^ 


= fT, 






f 


f 


f 


f 






= fE. 






f 


f 







(152) 



We turn to the terms fl55| [57j) . Applying the ansatz fll43p we derive 

JO.e.m. = \ F^" e,,^p [P^ (ec + eQQ)P-P(ec + eg Q) P") 

+ iF^^P^[-ec + ~eQQ)Pu + --- . (153) 

where we identify = My for convenience. Like in (11441) . we drop in (11531) 
additional irrelevant terms. The interaction terms displayed in (I153P are readi- 
ly reproduced using the heavy-quark multiplet field H . We follow [^8|50] and 
write 

Ce.m. = +^F^,trHQa^,H+^ F^,, tr a^, H H , (154) 

where the term proportional to eg respects the heavy-quark symmetry, but the 
term proportional to ec breaks it. We expect jegl ^ \ec\ since the parameter 
6(7 should vanish in the heavy-quark limit. The requirement that (I153p and 
(I154p agree yields the desired relations 



ec = ec , eg = eg . (155) 

Note that the empirical values (I140p confirm the expectation Icgl ^ |ec|. 

Finally we provide an analysis of the effective resonance interaction terms (17^ 
IHOl) . Assuming a decomposition analogous to (11430 for the resonance fields R 
and R^i,, we introduce a resonance multiplet field 



^=2(1 + /') (757^^^ + ^+)' 

^ = 70 5t 70 = (4,^ 75 7m + Ri)l{l+^), v>^K = 0, (156) 

where the field S transforms like the field H under a spin rotation (see I148p . 
At leading order in a heavy-quark mass expansion, the interaction terms (174t 
ISUil are described by 
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£ = ^ tr (5 7m {drV^n H +H^, {d^V^^^) S 

+ ^tr[Sj,^^{d>'^)H + H 75 7m (5^$) ^) , (157) 
which imphes the identifications 



Appendix C 

All tensor integrals will be decomposed into scalar objects of the form 



lab = J -(2^ Sail) S^il + p) , lab = -I J Sa{l) Sf^{l + p) , 

f d'^l 

Jabc = J Sai^) ^bi^ + Q) Sc{l + P) , 

Jabc = +l j ^ Sail) S,{1 + P) S^l + p) . (159) 

We introduce the two integrals Jabc and Jabc even though Jabc = Jcba as this 
redundancy will be used for a consistency check of the numerical simulation. 
The integrals lab and lab are ultraviolet divergent and read (see (fT5|l ) 

\2ml-ml 2 s J \ml) J 

2 s - {nia + nibf ( , ,2 



Vab = - (rua -rrih) ) , (160) 

where the logarithmic divergence sits in the subtraction term labi^)- The dif- 
ference Iab{s) — /afe(O) is finite. According to the discussion of Section 3.3, we 
define the renormalized expressions 



lab 



Iab{P^) - Lbif^ll) , Lb Iab{P^) - LbifJ'll) , (161) 



with the matching scale /im- The integrals Jabc and Jabc are finite. They may 
be evaluated in terms of their dispersion-integral representations 
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o'afec — / — 2 ~ ' '^"''c ~ / 2 ~ ' V-*-"^/ 



with the spectral densities 

^ log + yjplck^) - log (-/^lafec - V^L^ 



] ~Pac 4 "^fe 




2 



s-2p2, \Jml+pl^ + ^ml+pl^, (163) 



We wrote the logarithm so as to ensure a smooth spectral density. It is impor- 
tant to impose the dispersion-integral representation in terms of the proper 
variables, i.e. keeping = and = [p — qY fixed. This implies that the 
spectral densities have an implicit dependence on = — 2p-q. The alterna- 
tive representation of the integrals Jabc and Jabc derived through an application 
of Feynman's parametrization is also useful. We have 



J, 



abc 



Q[z2 - zl] e[^i - zl] e[zi - Z2]dzi dz2/{l6n' 



(m^ - zip"^) (1 - Zl) + ziml + 2{l- zi) z^ {p ■ q) + Z2 f4a 



r2\ 



J ^ _ f Q[z2- 4] Q[-gi - ^i] Q[l - ^1 - Z2] dzi dz2/{16n'j 

J {ml - Zl p2) (1 _ Zl) + ziml + 2zi Z2 {p ■ q) + ^2 /^L ' 

where we used = and = ml — m^. Depending on the values of the 
various mass parameters either (I162p or (11641) may be more economical in a 
numerical simulation, li p^ > (ma + m^)^ for instance, the representation (11621) 
is advantageous, the integrals being complex. We performed numerical checks 
using (I162p or (I164p and verify that the identity Jabc = Jcba is indeed satisfied. 



Appendix D 



We provide explicit expressions for the contributions of the (p D* and K*D* 
channels to the decay amplitude fl85l) . The latter are linear in the resonance 
coupling constant gn introduced in flHOj) . We obtain 



^ K-f A i- = YJ^{9T K%t (P' P) + ~9v {Kdi^P^ p) + (P, P)) } 
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+ 



+ 



+ 



9h 



-9v 



egH 
2/2 

9h 



{gr El%%, (p, p) + ~gy {E^^t* (P, P) + E-.%%. (p, p)) } 



2/2 M2 



{gTEl%D'ip^p) 

gy [E-J',UP,P)+E--%%*iP.P))} 



egH 



2P 



+ 
+ 
+ 



ge gu 



{gr E1%'*^, {p, p) + ~gy {Ett* (P, P) + Ef^U. {p, p)) } 
{cT E^.'S.K* (P, P) + e> {Ef4' (P, P) + ^TS'*^. (p, p)) } 
{ESe^.fep) + 2ESe^^fep)} 



[gTFf-,.{p.p)+gEFf-,.ip^p)} 



|ff^{F;5e),v(P,P)+Fi5e),v(P,P)} 

9p9H r pa/3,/. . X pa/3,/. . 
36 /3 my ^ ^ IP' PJ / ' 



(165) 



with 



+g-.. {^f (p + ^r'^ip + + 5f '"'(p + (p + 0}} , 

<i(p, P) = -'J(^. ^'^-^(O {p" '"'(p + + ^7^" 5f '''(p + /) 

+r {^f ''(p + sr'ip + + '""(p + srip + o}} 
<i(p, = - ^ / ''(0 5r^^(p- + /) 

+p-sr'iP + i) + g'''S^%p + i) 

+r {^f(p + 1) srip + + ^r"'(p + ^r(p + o}} , 
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{stiP + sH'ip + /) - s^^ip + 1) sr{p + /)}| , 

|r {Sf'^ip + /) S^'ip + /) - Sf^^ip + I) S^'ip + /)}| , 

|r g« {srip + 1) srip + - s^ip + /) + /)}| , 

x{p + lY{p + irS,{p + l)Sl^{p + l), 

x{p + ly ip + I), p^ S,{p + /) S^^^ip + /) . (166) 

The evaluation of the decay constant d^+^^i- defined in f l87p is done by con- 
tracting the gauge-invariant tensors [HH ES]) and fll66p with the anti- 
symmetric tensor 



We derive the required contractions in terms of the master loop integrals 
lab, lab, Jabc and Jabc iutroduced in (11591) . Following the arguments of Section 
3.3, reduced tadpole integrals are dropped systematically. Using the notation 
p^ = Ml and {p — qY = Mj, we have 



8(p-g)M/P;;, 
-Mj(2 



2 d(i-) 



/3,M 
2 



a6 



mi. 



2 m? + Mf 



2 Mj) 



Ml 



'ofe 
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+ 2 Mj (ml -ml + M}) Mf U + Ami Mf [Mf - M}) J^ab 



48 {p ■ q) Mt P%-1 AT = - [m] - Ml 

+ 3Mf + {ml + mt) M} - 2 (m^ - 



3 [ml + mt) - Mf) Mf 



mi 



l-ab : 



(Mj - 2 Mf Mj - Mf 



m„ 



+ Mf (m) + 2 Mf Mj - Mf) + (-M) + 2 Mf M| + M;^ 
- 2 M| (-m^ + + mJ) Mf 7,6 + 4 M| Mf (m| - Mf ) J^^b ^ 



'aft 



8 (p ■ g) Mf Pi;-) C, 
-2 



/3,^ ^ abc 



2 (m„ -m^ 



MJ 



GMfMf 



3ml - 2ml - ml + Mf] M? Mf - 4 M? 



2 MjMf 
AMjMf 



2 ml- + ml- Mf) Mf + ( 
2M;+(m^ 

(1-) ^af3,fi 



m. 



m, 



mt -ml + Mf) Mf 
)^ Mf + 2ml Mf Mj - Mf 



12 (p ■ g) Mf P(;3-) C/^/ = - [ (4 - 2 (mJ - M,M, 



Mj - Mf) (m| + 3Mf) - 3m^ (m) - 3 Mf Mf) ) Mf ml 



(ml - Mf) (Mf - Mf) (-3 Mf + MfMf + 2 ml Mf 



+ 3mlMf 



Mf 



Mf) (2 Mf 



MfMf 



mlMf 



QmlMfMf 



-3miMf 



Mf 



-2 ml + 3 ml- ml + Mf 



m. 



QmiMf 



mi + M M. 



m, 



'Mj+{ml 



I be 



mi 



Mi 



2 ml Mf Mf 



miMf 



J, 



abc ■ 



12 (p-g)Mfpf-) 
-3(Mf- 
{ml 



4 D(,l-i pal3,ti 
Q-P,^ abc 



6 



mt) (ml 



M?) Mf 



Mf 



2 m: 



m„ 



m. 



hmi + 3mi]M, 



Mt 



+ 2 (Mf - Mf) (ml - 2 (ml + Mf) + - 2Mf + m^Mf 
- (Mf - Mj M,)' (2 ml - (3ml + Aml + 10 Mf 
+ 2(m2-Mf)' + 3m2 (m^ + Mf ^ 



m„ 



3 ( Mf - Mj) Mf 



3 (m^ + Mf) Mf [mJ (-m^ + 3 - 2 + Mf) 



— —m„ 



mi + Mf) M^ 



Lab 
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6(ml + Mfj Mf [Wjml - {2mlMj + (mJ - 



J, 



ahc ■ 



9QM}Mt{p-q)P%-}E:^^^ 



-3 Mj 



A {{ml 



m; 



-ml + ml + Mf) {Mf - MJ) {-ml -3ml + 7Mf) 



+ 



+ (Mj-Mf) {-3M^-2{ml + 3ml)M^ + {ml- 
- {5Mj-{7ml- 



'a6 



M«- 17 M^m^Mf 



M| + 2 {ml - rr 

- Mj) ml + 5ml + 5 Mf - 22ml ^f) ^1 



Mj{{2 



19 Ml 



- 19 + 2 + 23 mt Mfj M^ 
6 ml M] Mf (m| - M^) [ - M/ + (-5 ml + ml + Mfj Mf 



lab 



M; (9m 



3 ml 



AM] 



J, 



abb • 



96 M? ■ q) Pi^J E^J^ = 6 M, 



4 ml 



Mf 



^{ml 



2mi ] M? 



5ml -Ami M^^^ 



m„ 



+ Mj {2ml + {Sml - 3 M 
+ {ml - ml) Mf {-ml - 5ml + M] 
+ 6MjM^ {-ml + ml + Mf) (2 ml + 10 

'2m2 + 10m? + M?-3M2 



10m^ + M;-7m^M;jM/ 



ml + M}-3Mf 



UmlMjMf 



Mj 



M-f]J, 



abb ■ 



96 M| {p . q) P^^,-} <i - o 



3M? 



+ (Mj-Mf) (,„4_^4^^4_2 
+ {Ml - Mf) {ml + 3Q mlml + 11 ml + M^ 
-2(ml + lQmt)Ml)Mf 



2mt] M- 



+ 2 ml {m] - Mff {-5ml + A mlml + ml + 9Mf 



2(5ml + 2ml)Mf 



Lab 



+ Mf Mf [2 Mj - {ml + 25 m^) Mf + (-m^ + 56 ml ml + 17 m^) Mf 
-5 ml+ {a Mf -Sml^ ml + 13 ml + Mf -2 ml Mf) Mf Mf 



12 ml 



-5 ml + 4 ml ml + m 



Mf 



{Mf -{5ml + 11 ml) Mf + A {ml 



mr 



M 



lab 



+ 6mlMf Mf {Mf - Mf) [a {5ml + ml - Mf) ml -2 {Mf - Mf 
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J, 



abb ■ 



48 Mf Mi P'T. E 



Mi 



f al3,f^ ab f 

- {3 ml + 2 {12 ml + M 

- 2 Mj (m^ - 9 ml ml - 2ml 



3M,'° + 2 (3m^ + 6m^ + M}) Mf 



m:,+9mf-Mt + 5mtM 



2 I 2 



M] 



+ M' (7 ml ml 



4 2 

'^b 



m, 



+ 2M^ml 



':+{5ml 
(^ml + m 



Mt 

6 mt ml + 5 mi 



2 (^m^ — 77 

4 (m^ - Mj) - 3 ((m^ - 2MJ) MJ + 2 (m^ + mJ) 
+ Mj (6 ml + 25 Mf - 13 Mf) + (ml + Mfj ( - 2 M' 
+ 2M^ Mj + ml (mJ + 2 M^) +ml(AMj-7 Mj 
+ 6 (3 + - M^) (m^ - MJ) 7,^^ , 



/ 

lab 



Mj]M, 



m„ 



2AMJ P^^-^E""^'^ 



ab 



-Mj 



3 ml 



mi+mi + Mf]M: 



(Mf - M|) (m^ + (m^ - 2 M^) - 2m: 



3 H 



m. 



Ml 



m^ + (m^ 



M 



Mf + miMi 



M 



-3Mj[ml-{ml^....f 
+ 6 ml Mj Mj (Mj - M]) J^bb 



2M?lm? + M5-3m? 



Mt 



lab 



lab 



2mn 



M' 



— 6 [m 

+ 3 (ml- ml) + 3 (mJ - 4 m^) MJ 
(-5 ml + A ml ml + ml - 8 Mj + (l3m^ 



'ab 



m, 



5M^^-(7m^ + m2 Mj + 2(m2-mn 



lab 



m. 



2M/ + M/) J,,, 



-6mgMjMf(Mj-Mf) (m^ + 
48 Pi;,-) 4f ''^ = {Mf - M|) [3 Mf - (6 [ml + 



mi 



M? 



+ (3 (m^ + 10 ml ml + m^) - 4 (m^ + m^) Mj) 



5 m^ - 6 m^ m^ + 5 ml) M} M^ 



-2[ml-ml 



12 MfP%-}Ef^^ 



-"ab 



M| (m, 



— 5 m„ + 4 ( m 



Mt 



m„ 



m. 



'aft : 



^^-^i ^aP,fi ^+,ab 



M] (M^ 



m^ + (4 



m. 



2M; 



m„ 
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5mt + M^ + AmiMf 



24 p(i~) 



-Mf (m= 



5 + 4 (^m^ 



lab 



2 Mf - 3 Mj 



2m2 + 3 + M]) + [ml 



m. 



M 



+ ml Mf [(3 ml (3Mj - ) Mj + 4 (m^ + Mj) (mJ - M/ 



rrtu 



2{M]-Mf 



6m^M;-2(m2-Mj) (m| - M^ 



+ 6 Mj Mf 



+ 3 m; 

,2 _2 /I ^2 
r2 



f-ml 



M- 



M'jml 



2m^Mj 



l-ab 



M] 



M 



m. 



+ mZ M: 



J, 



abc ■ 



12 Mf ■ g) P^^,^ . 



3m^M^ 



2 Mf - 3 M| Mf 



M^ 



M} 



-m„ 



2ml + ?>ml + Mf] M, 



M}Mf 



2{Mf -M- 



3 (3M| - Mf) Mj + 4 (m^ + M|) (m| - M: 

2(m^-M|)'(Mj-M, 



ml + Gm^^M} 



lab 



+ 3 M| (mJ - (Mf - Mf] 
6mlMjMf[Mfml- (2 Mj + (mJ - Mf 



+ mlM 



m 



Appendix E 



We provide the contributions of the (t)D* and K*D* channels to the decay 
amplitude M^+^f^Q. introduced in fl96|) . All terms are linear in the coupling 
constant gu (see flHOj) ). We obtain 



- 8 j2 \^-,<A02 



3 J2 V+,'f>Dl 

9h 



16/2M2 



[ec-e+ 3 J 
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+ 



+ 



9h 



{qe F'^f.j,, ip, p) + ~gr F'j^i^. {p, p) } 



(2ec + eQ/3) qh 



9v H^^DD*{P,P) 



4/2 M2 



+ r^:TrTT^^ Qv HT^j^ij. [p, p) 



eA 9p Oh 



ca 9p 9h 



24 /3 my 



18/3 my 



(169) 



with 



+ p'p- g-^^ {Sf'^p + /) Sriv + + Sf^^%p + /) Srip + /)} 
E'Jt{p.p) = +^ J -^,e^',fe.-,.fSf'^^l) \p'p-Sl'^^{p + l) 

+r g'"" sr{p + /) + /^p" sr{p p^p^ sf'^p + 1) 

+ p'p'^ g,. {Srip + /) Snp + /) + Sp'-'^ip + /) S^ip + /)}| , 

n::'ip^p) = +^ J -^,^^-,f^.-p.fS7{i) ^p'g'-sf'^'^{p+i) 

+ sf^^%p + i)st^~^{p + i)]Y 

{^r(p + sr^'^p + + + sr^p + o} 
p>'' q. {sf'^'ip + 1) sr\p +1)- sf'^ip + /) sr^P + 0} 
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{srip + 1) sr^ip + /) - si^ip + /) sr^P + /)}} 
x(/+p)'pp5fe(p+0'5r(p+/), 

X (/ + S,ip + /) + /) . (170) 

The decay constant di+^^o- defined in ( !97|) is computed by contracting the 
gauge-invariant tensors (pHI - llOll I104p and (11701) with the antisymmetric tensor 



p(o-) Iff P^q'^X r p^q\ \ 

= "2 [{9f^- ~J7^\P^~ W " TTVi^^j ■ ^^^^^ 

We provide the results in terms of the master loop integrals lab, lab, Jabc and 
Jabc introduced in (11591) . According to the arguments of Section 3.3 reduced 
tadpole integrals are dropped. Using = Mf and (p — qY = M?, we derive 



16(p-g)Mjp(°-jAX" 



2Mj 



+ 



Mf - ( 



— m 

2 



Lab 



Mf -2(ml + ml) + (ml - mlf Mf 



+ 2[{ml- 
2 



m. 



7 



+ AmlM}Mf 



m„ — mi. 



m„ 



f 



m. 



(Ml 



-M) 



J. 



Mf 



Mf 



l-ab 



aab 



(Ml 



(Ml 

Mf + (ml 



3M? 



m, 



Mf 



Ml 

, 2 



m^ 



lab , 



16 (p ■ q) Mj PfJ BT' = M 
- (Mf - (5ml + 3r 



m„ — m. 



M) 



8ml + Mf 



Mf 



2 mt 



mi 



Mf 



Mf -2(ml + ml) Mf + (ml 



mi 
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(m; + {ml - mt) MJ - 2 {ml 



-2Mj 



m„ 



hmi]M 



'/ 



2mlM]MHM^,-Mf] (-2 m: 



f 



m„ 



m, 



M 



Lab 



2ml-M^j+3M^)jatb. 



1 9 P^*^^) Uf^'<^f^ 



Mt 



M] 



m„ 



2ml + Mf + 7mlMt 



lab ■ 



2 M?) ml 



3M! + 3{ml-miyM^ 



+ {Mf-M]) {Mt+{ml + ml)M^-2 



m„ 



mi 



lab 



AmlM] 



2 Mf + 3 M; Mf 



16 {p ■ q) Mj P'^'J C,,, 

+ Mj (-3 ml + 2ml + ^fj + {ml 

-4mlMjMf[{2ml-3ml + ml-M])Mf 



m. 



'be 



%mlM)M} 



mlMj - {ml - mlf 



2ml Mi M} 



J, 



abc ■ 



12 {p-q)M]M^Pl^'Jc!:if 

2 



2iMJ-MfM- 



m 



m„ —2m 

Mj - M^)' (mJ + 3M2) - 3m2 (m) - 3 M] Mf) ) M^m^ 

3Mf + MjMf + 2m2Mj) 



7 



6m^M2Mf -3m2Mf {m^ - Mj) {2 Mf - Mj Mf - ml 



3mlM}Mf 
QmlMfM^ 



{2 ml - 3 ml + ml- Mfj Mj 



m,. 



I be 



mlMj 



Mi -2mlMt Mi + mlMf 



J, 



abc ■ 



48 M| P^^J Di:ir^ = 2 M| [ - 3 (-m^ + + M^) 



+ 



6 (m^ - 2 {ml + M 
- 6 Mj (m^ - (m^ + 2 M 
+ 12 ml M] {Mf - M}) Jabb 



2Mt]mi-2ml + Mf + miMt 

, 2 



lab 



m. 



ml + M) 



Mf 



3miMj]M, 



lab 



(172) 



,(0-) 



^+,ab 



2Aml {p-q)M^P^ 

+ 3 (Mf - Mj) {ml + 20 ml m, 



24 m, 



m„ 



11 mf 



-mf) -Mf]M\ 
-Mf 
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2 {ml + 2m^) uf) - 2 (mJ - M^)' (m^ + (7m^ - 2M^) 



Lab 



ml + 3 ml) MjMf + 4(m^ 



2 Mj M?^ 



-A'/) + 10 + (ml -8 ml ml -9 mt) Adf) Mf 

+ M} (-Mj + (m^ + mt) Mj + 2ml (ml - ml) Mj) 
+ 4mlMj (ml - 2 (ml + Mf) ml + {ml - uf) ' 



+ (jnl — n 
+ 6mlMfMf [m^ -M 

^', + Mj)Mf+ml{8 



8m^ + 7Mf m^ 



ml + 3M'j 



b 

3M? 



2 + M/ 



J, 



abb 



2Aml{p-q)M}P^^jE>!:ll = M^ 



24 mf 



ml -ml 



-Mf] M 



3{Mf -M] 



mt + 20mimt + llmf + Mf 



2 [ml + 2m{) Mf) + 2 (mJ - Mf)' (m\ 



m. 



2Mt 



2ml + Mf + 7mlMt 



Lab 



m„ 



3 ml) M]Mf + A (m; 



ml m.b 



-M' 



/ 



10 ml Mf 



m„ 



1 m, m„ 



^ (2 Mj Mf - m;^ 



+ M;(-M? + (mf + mf)M^ + 2m: 



m„ 



m, 



Mj]M, 



+ AmlMf (m^ - 2 (m^ + Mj) ml + (ml - M' 



m„ 



m. 



M? 



lab 



+ 6 m^ Mj Mf (Mf - M 



(m^ + Mj) Mf + m^ (s 



im. 



m. 



-7M]ml 
3 M j - 3 M; 



2Mj 



M/ 



J, 



abb ; 



m„ 



1 9 M^ P^°~^ pM-"/3 



+ [2 [m'^ 
+ 2 (ml 



m. 



m^ 



(Mf 



Mj 



M] 



TT^a 
lab ■ 



mi 



7Mf + 
Mfj Mf 



b(ml + 



mi 



-Mj 



Mf 



(173) 



24 {p . g) M| P(%) F^;,t = 2 M,^ 



+ 3mf 
+ K 



7 ml + 9 ml 

-Mr 



12 m; 



ml + ml + Mt] M 



7M,- 



7m„ 



(5 



(Mf - M 
11 M^ 



M/ 



m 



m„ 
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- 6 [ - 2 + 4 (m^ 

+ [Aml-Aml + Mf) (-m^ + 

+ ( - 3M) + (m^ - 9 mfj Mf + 2 (ml - m 
■ I2mlml M] Mf (Ami - Ami + M] - b Mf 



Mf) Mj 
Mf 



Mf 



J. 



abb 



24 ml (p ■ q) M] PjfJ F^_%^ = 6 ml M} [a (-ml + 
+ (Mf-Mj) (7ml + 9ml~7Mf)Mf 
+ Mf)(Mf-M] ^ 



m. 



M?) M? 



+ [m 



ml 



6 ml [2 (ml - 2 (m^ + uf) ml + (ml - Mfj')Mj 



- [Ami - Ami + M]) {-^a + ^ 



Mf) Mf Mj 



3Mf 



il-9ml)M]-2(ml 
+ l2mlM]Mf (M'j-Mf) (-4m2 +4r 



— m\ 



Mf 



l-ab 



■Mf + 5 Mf 



J, 



abb 



(174) 



2 p!.'J K::^ = ml (-2 Mf + Mj Mf + yn^ - ml ) MJ ) l^b 
+ ml (-ml + ml + Mf) Mf U + 2mlml Mf (Mf - Mf) Jabb 



ml - ml 



M] 



Q Mf P^'J Ft!:-J = 



3 Mf - (a m 



9ml + AMf 



+ Mj(2ml + 5ml + Mf) Mf-2 (ml 
+ Mj(2ml-mlml-mt + 
- 3 ml (-ml + + Mf) Mf J, 



2 I 2 



+ 



6mtMf(Mf-MfJM 



Mf 

'm} 

Mf] I, 



lab 



(175) 



12 {p . q) M] P^^J GXf 
+ Mf(-ml-2ml 



2Mf 



?>mlM} 



3ml + Mj] Mf + [m 



3M]Mf 



m. 



Mj 



+ Mf [(3 ml (3Mj - Mf) Mf + A (ml + MJ) (mJ - Mf 



2(M'f-Mf) mt-6mtMi -2 



+ 3mlM] 
GmlMjMf 



Mf - m 
M'jml 



f 

Mf 



I) {Mf - 
2mlMf + IM 



lab 
2 



U) 



m. 



Mt 



Mf 



m. 



m^Mf 



J, 



abc 
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48 {p ■ q) M} p(°-j H^X/ = -2 Mj [6 {ml - m^) {ml + ml - ) Mf 
+ 3 (Mf - M]) {ml +{-3 ml + 3 ml -2 



{ml - Mf) {ml + Ml 



m„ 



Mf 



m: 



{Ami 



- 2 + 5 M/ 



2mt + Mf + 7mlM: 



-6M: 



{2Mf-3 MfM} + Mf) ml -2 {ml {3 Mf - 3 Mf Mf + Mf) 



Mf{2Mf-Mf) {ml 



Mf 



Mt 



m„ 



m^, 



Mf) {2Mf 



-3 m, 



2ml + 3 Mf) Mf Mf + {Mf - mf) M^ 



l-ab 



- 12 Mf Mf 



: {Mf - M 
+ ml {ml + ml-Mf 
— mA M^ Im 



m,. 



Mf 



mu 



Mt 



2 



2 ml- M? 



Mf -Mf 



Mt 



Jabc • 



(176) 



Appendix F 



We derive tlie contribution of the (pD* and K*D* channels to the l"*" ^ 7O"'" 
process. According to (11071) it suffices to evaluate the contractions of the four 
loop tensors introduced in (11081 11091 IllOp with the projector 

We have 



2Mf 



P^^ A^t'X = 2 ml Mf 4. + 2 ml {ml - mf) Mf J, 
2 ml + {Mf - Mf)) Mf + {ml - ml) {Mf - M, 



{p .q)Mf' ^^-^ 
- '{2rnl 

24 



abc 



lab 



f p(0+) _ ,^2 

^u.a0 ^+.ab — f 



-2Mf+{ml-llml-Mf)Mt 



m„ 



5 (8m^ 



-2f2m^ 



Mf)ml + 7{mt-mlMf))Mf 



t-mimi-mt)Mf 



lab 



+ 



3 {ml -mt + Mj-2 {ml -2ml) Mf) M, 



-6 ml Mf (7 ml + ml- Mf) Mf\ hb 
+ 6 ml Mf Mf (3 ml + ml- Mf - 2 Mf) {Mf - Mf) Jabb 
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^^f p(0+) ^ti,al3 _ J.. 2 



+ [ [ml + M}) (Ml - Mf) - ml (? Mj + M, 



+ 8 Ml (mJ - ) J.^, , 



9/1 p(0+) nMiO/^ _ 



3m^ (ml + 3m^ + Mf) (m^ - uf) Mf 



- (mJ - M^y (ml + (Ami - 2Mf) - 5m^ + M/ + 4m^Mf) 



- 6 



(ml - 



mi]Mf 



+ 3mlMl 
+ 2m 
+ 6m^Mf 



2mlM'j 



3 ml 
he 



Mj) (mJ - Mf) - 



m,, 



ml+ml + M^ 



M"' 



J, 



abc 



(178) 



where we use = Mf and (p — Q')^ = Mj. 



Appendix G 



We collect explicit expressions for the contributions to the decay amplitude 
(11121) which are linear in coupling constant qh (see (l80l) ). We have 



+ {F^TDf{P,p) - EfSf/.{p,p) + F^%%^,{p,p)} 

+ Sf^M^ [ec - e + ^] {gr F^'^^^f{p, p) 

+ ~9v {Ff^r'{P.v) + F^_%f{p,p))} 

(^|er^(p,p) + ^-"it'Se(p,p))} 
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9h 



[cT E^'^l (p, p) + e> {Ef£:t (p, p) + Ef4f^{p, p)) } 

'V 

- 1^1^ [GtrAp.p) + GtrAp,p)} 

eA~9P9H {Gf,i;f^p,p) + p)} , (179) 



64 /3 my 



+ 9% P" St^'ip + l)+ g% g^^ Stip + /) + 9% p" ^^^'^Ip + 
+ 9% 9-n. P'^ {Stip + /) Sr'ip + + St'-^ip + I) S^%p + /) } I , 

+ r + /) + r sr'^'ip + + srip + o 
+p>" {sf%p + /) 5r'^(p + /) + s^'ip + /) 5r(p + /)}| , 

E^-ir'ip,p) = +^ J (0ie.;3p^-^r'^(o{^?">"^f'''(p+o 
+ f g^^ sf\p +i)+p^ sf^'-'ip + /) I , 

2^ J ^,^^f,fp''S!'^^l)p^Sr{p + l)qr, 

^ j^.^.-,fp^sf^''^^)p^snp+i) 

E^-T'iP,P) = ^ I 7^,^^^fp"Sr^{l)p-sf'^{p + l) 



ET'''''iP,p) = 

Kt'^'-^p.p) = 

Eiir'ip,p) = 
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Ft'^'-'iP.p) = / |p>" s^^,M sr'iP + /) 

+ 9% 9-.. P° [Stip + /) S^^^p + l) + S!'^%p + /) St^'^ip + /) } 

+p>" {sf{p + 5r'"'^(p + /) + si'^%p + /) + /)} 

g% q^p- [stiP + /) Sr'^{p + 1)- St^'ip + /) sr^P + 0}} 

r g« {^f '^(p + sr\p +1)- sf'^ip + 1) srhp + } J 

^+,ab \P^P) - ^ J (27r)4 a^P ^'^'^ 

I rp'^ g« {5r(p + sr^p + /) - 5r(p + o sr^p + o}} 

x{p + ly ip + z)'^ Skip + 1) s:p{p + /) , 

ap+ ly ip + s,{p + /) s^^ip + /) , 

X (p + + If StiP + /) S^'^'^ip + /) , 
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■p e. 



a0 



(180) 



According to flllSp it is enough to compute the contractions of the tensor 
integrals (11801) with two tensors 



p(2) 



^ P" <f {^ara^, Pp Qa P/3 - e^ra^, Pp 9/3 Pa 

- ear/3M P^ Pp + e^^p^ pa 9/3 Pa } , 
~\P'' 'f {ear^a Pp P/3 " ear/x/S Qa Pp Pa 

- (^ar^ia Qp Pa Pfi + Qp Pa Pa] ■ 



(181) 



We estabhsh the results 



af),^,ap +,0.0 



al3,fj,,al3 



12 ( 



m? - ml + M^) (ml + M^) 



(p-q) 

+ Q{Mf- Mfj (ml + 2M^ml + 5M^ + ml(Mf -1 
- (Mf - M^f (19 Mf - 8 (ml -2ml) Mf + (ml - 
-Q(ml-ml + Mj) (ml + Mj) (m] + M^) 
- 12 ml (ml + Mf) (mJ - M^) (mJ + M^) J,,, , 



48 -^pW 4 

afJ,^,afj —,ao 



al3,fj,,al3 



(p-q) 



12 (M 



2 



6(M^-Mf) (-mt-AM^ml + hMt + ml (ml + M^)) 
- (Mf-M^f(-17M^ + 4(ml + ml)M^ + (ml- 



ml 



Lab 



-6 (m 



M] 



-m„ 



Lab 



12 (m^ - Mf) Mf (mJ - Mf) (mJ + M^) J^.^ 



24 ^P^^ 



.(1) 



(p-q) 



_ _ Aa|3,^J.,af3 
3 al3,fi,ap ^+,ab 



Mf 



2 (ml- 2 (ml + 2Mj 



(m^ + (4m2 - 3 Mj) - 5m^ + M|) M/ 
+ 2 (3 - 7 + 4 m^) MJ M^ -(ml-m^f M) 



lab 
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im? M^M' 



M2±mg (p-g) 



TjO!P,fi,aP 
aP,fi,al3 ^±,ab 



M? 



24 m^M^ 



m^ +m^ M/ 
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(m^ - Mj) {ml + 10 ml - 11 + - 2 (m^ + ml) uf) 



M' 



Mj - Ml 



19mf M 



2 (m^ - 5 mfj Mf + (m 



m„ + 10 m^, m„ 

-4 

+ + (m^ + 7m2) M} - 2 (m^ - 



11 



+ 2Mj 



3 ( -m^ + Qmlml-bml + M^f -I2ml 



QmlM] Mf (m} - M\ 
+ 3 [m] - M, 



2Mj + MfM'} + Mf 



mi 



[m^ + imI 



J, 



abb ■ 



MfMj 



48 



M]Mf 



''J_ p(l_) r)Oil3,tj.,al3 



10 



2 — Ami 



m„ 



(4 



M: 



5 



M) 



Mf\ Mf 



10 ml Mf) Mf 



mr m 



l + {3Ml 

- ^b) ^/ - 



mi + 5Mf]]M, 



3Mj 



m„ 



hmlMfj Mf 



m„ 



mr 



Mf 



lab 



+ {ml + Mj) Mf 



2Mf 



ml - 17 ml) Mf + {ml - mlf Mf 



Mf 



7ml)M^f + 2(ml-m^) W, 



+ 6m^ Mf {ml + Mj) {Mf - Mf) [- ml + ml + 2 Mf - M, 



J. 



abb ■ 



12 M? 



>(i) 



(p-q) 



B 



3 oi/3,fi,al3 ab 



a0,fi,af3 



Mf 



2Mf-{ml + 



{ml-miy 



Mi 



12 {ml-mlmbf Mf 
12 ml {ml - ml) {-ml + ml - Mf) {Mf - Mf) Mf 



192 M? Mf P^^^ „ CI 



+ m2(Mj-Mf) (-5Mf + 4(m2 + mgMf + (m2 
+ {Mf - MfY {SMf - {7ml-5ml + 5 ml) Mf 
+ (5 + 8 + - 2 (3 + 8 ml) ml) Mf 

- {ml-mlf {ml + 



mi. — m. 



Mf 



+ {Mf - MfY { -2Mf + {5ml + 9ml) Mf 
-2 {2ml - ml + 3m^) Mf + {ml - ml 



<Mf 



2[Mf 



8MfMf-AMfMf + Mf 



mi. 



- {Mf - Mf) (8 Mf + 5mlMf -3{3ml + 7 Mf) MfM 



lab : 



Mf 
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+ (Aml + Mf) Mf) ml-Qml Mj (m} + M] 



+ [-ml + ml + 5 Mf) Mf) + ( 
-2M'fml + bM] [Mf -2M] 



3mt + mt + Mf] Mj - 2 (Smt + mt + 3 Mf] Mf M 



M? 



M'jml 



7 



Lbc 



+ 6 ml Mf Mf [2 [ml - mlf M^ - 2 ml (mJ - Mf 
+ [ml - ml 
- [Mf - M 



m„ 



m„ 



my 



M] 



m„ 



3ml + 2ml + ml-3Mf) [Mf - Mf 
+ 2ml-2Mf)ml + ml [Mf - 3 
Mf 



m 



M? 

2 



-Mf 
mt. 



o 2 2 

3 m, — m. 



J, 



abc ■ 



192 M^fMfP^'^ ^C'^'^'tf 

J * ap,^,afJ —,aoc 



Mf 



12 m 



m 



Mf 



- 12 ml [ml - ml) [-ml + ml 
+ ml [m 



3Mf 



Mf - Mf] Mf 



m„ 



2 _ 



Mt 



7Mf + 8 5m 



3 ml 



Mf 



— m„ — 



+ [rnl- 



4 ml ml — 3 m\ + 2 [mi — 2 ml 



+ [m„ -m^ 



Mt-[Mf-Mf) [Mf-[ml-lml + 3 

^9 Mf 

ml - ml) ) Mf 
Mf - MfY [-2M^ 



m. 



Mf 



[^l 
Mf 



m„ 



m. 



+ m 



m. 



Mt 



+ mtM, 



8 MfMf -AMfMf + Mf 



-2[m\ 

-Qm^Mj [Mf + Mf) + [Mf - Mf) [a Mf 
+ ml [-5Mf + 9Mf Mf - AM, 
+ [ml - Mf) (m 



m. 



33 Mf Mf - Mf Mf 



m: 



[2Mf + 5 Mf Mf 



m. 



[Mf-2Mf)) 



-3mlMf 



-3 Mj -QMfM^f + Mf 



m,. 



+ [Mf - Mf) [7Mj + 5 Mf Mf + mf [Mf - Mf 
+ QmlMfMt[2[ml-mtfMt + 2ml[Mf-Ml ' 



m„ 



my 



-3ml + 2ml + ml - 7 Mf 



I- be 



Mf - Mf) Mf 



+ [ml - ml) [Mf - Mff [AMf + [Uml - 5ml - A 



m. 



Mf 



my 



[^l - 



m. 



Mf-Mf)\Aml-3 



mi, m„ 
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J, 



abc 



192 Mf P. 



(1) 



al3,fj.,al3 ^abc 

+ 12 (ml - ml) (ml + 



-Mf 



12 [ml 



mi 



Mf 



m. 



m. 



2 Mj) M/ 



Mi - Mf 



MfMi 



Ml 

2 



2(ml-bml)M^ + (ml-m^ 



(m^ + (4 



m^ 



2M; 



m„ 



+ Mf -2(2Mf + MfMj- 11 M/Mj + 2 Mf 



m; - M- 



2MJ 



9MfM] 



11 M/M? 



+ (m) - 21 A'/J + 8M^))ml + 6ml Mf (mJ + M: 
il - Mf) (m] - Mff (2 Mf -MfMf + ml [Mf + 4 



m,. 



3mlMf 



-M^ + QMfMi + ^Mf 



m. 



,^ -r u ivii ivif 
+ [m] - Mf) [Mf -MfMf+ ml [m] + 3 



'6c 



+ 6MjM^[ (m^ - mt) (Sml + 2ml - 3ml + Mf) (m^ - M|) Mf 



+ m2 (ml + m2 + Mf) (Mf - Mj)' - 2 (m^ - m; 



Mf 



m„ 



mi 



Mf-M- 



ml + 3 ml 



Mf 



+ (ml + ml) (ml - ml) ) 



J, 



abc 



192 



P. 



(1) 



,D 



+ (ml - m2) Mj 



12 ml Mf 



Mf 



ml — 2 ml + ml 



Mj]M, 



+ 4 



3 ml Mf (Mf + Mf ) - 2 (ml + M|) (Mf - Ml 
Mf - M2) ' + 6 ml Mf - (ml - Mf) ' (Mf - Mt 



+ 3 ml Mf (Mf - ml) (Mf - Mf 



lab 



- 24 ml Mf [ ( 

+ [rr 
M'iM} 



mlMf 



m. 



mlMf 



f-Mf]ml 



96 



Mj 



Mf 



Mf 



17 mi -25 mi 



J abc 

15 Mf) Mf 



2 (7 M) - 3 (m^ + 7m^) M| + 4 (2 + 5 - 7 



m;^ 



+ (m^ + 7mi) Mf [3mt-3mi + Mf) Mf 

, 2 



13 Af* 



m; 



86 



- 2 (-5 M} + 4 {ml- 2ml) + {ml - mlf^ Mf 

+ 3M] {3ml + 2mlml-5ml + Mj -A{ml + 3 m^) Mj) Mf 
- Mf {-5ml - Uml ml + 19 m^ + Mf + 4 {ml + ml) Nlfj Mf 
+ 6m^ Mj (mJ - Ml)[ - AMf + {ml + 7 ml + uf) 
+ Mj (3 + 5 - Mj 



(p ■ q) aP,^,af3 + 



,af3,fi,al3 



-Mj 



J abb 1 

12 



Mrf-mi)Mt 



-3(Mf 



M'j-M'f 



5Mt + A{ml + 3ml) + {ml - ml 
AMf + {iml + mt) + (m^ - 



ml 



iab 



+ 



- 2 { - 5 + 4 {ml - 2ml) M^j + {ml - ^1 



Mf 



+ 



+ M^f{M'f-^{ml + ml)M} + 7 {ml-mt) ) Mf 

M^j {7 ml + 10 ml ml - 17 + - 8 {ml + ml) M}) Mf 
Mj Mf {m} -Mf)\-AMt+ {ml - ml + Mj) Mf 



l-ab 



+ 6mlMfMf {m^-M^ 
+ Mf{3ml 



5 ml 



Ms 



-AMt+{ml 

J abb 1 



ml Mf (I) 
96 -7^ — ^ P^^> 



(p-q) 



paf3,fj,,af3 
al3,^i,al3 ^-,ab 



-2 M? Mf 



- 12ml 



-ml + ml + Mf)Mf 



Mj - Mf) (m/ -2 {ml- 11 ml) Mf + {mf - ml) 



l-ab 



+ 



-M; 



Mf 



2 (ml -2{ml + MJ) ml + {ml 

-2Mf(ml-2 {ml + MJ) ml + {ml - M 

- 24 ml Mj {-ml + ml + Mf) M] 
+ 48 ml Mj Mf {Mj - Mf) J, 



Mf 



3 al3,fj.,aP ab 



(p-q) 



J abb 



— ml — 4 + 5 ml — 



+ 2 {ml- 2ml) M: 



l-ab 



3Ml 
(p-q) 



3 afi,ii,af3 ab } 



-2Mf 



m„ 



ml) Mf 



{ml -ml)' 



l-ab 
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iaf3,fi,af3 



M] 



5 m„ 



18 ml ml + m'D + (ml — m 



-3Mf + (7ml + 5 

lab 



m\ 



Mf 



96 



E 



, ab 



-Mi 



2 I 2 

m^ + m^ 



m: 



al3,fi,aP ab 



(p-q) 
- (2 ml 



M) 



lab 
2 



{Mf-Mfy{ml 



3 Mf) ml + {mt + 9 Mf ml) ml + 2 (Mf - ml M, 
Mf - Mfj (ml + 18 ml m^ - (l5 + 4 + 3 M 



-2 m 



2Mfmt 



IMfml 



+ 12 ml (3 ml + ml- Mf) (-ml + ml + Mf) M, 
+ Mf [ (ml -ml + Mf) (m^ + (l9 ml + M 

- 2 Mfj Mf + (2 - 4 (m^ + 10 + m^) M 



Mf] M 



ml + Ami - 2Mj 



2 m; 



2 

Mj — 3 (ml — 6 m^ m^ + m^ m^ + 4 m| 



+ 2 (m|| - 3 m^ m^ + 2 m^ - (2 m^ - 5 m^) 



+ (m 



12 m^ ml 



m, 



Mf M. 



lab 
2 



+ 6 m^ Mj Mf (Mf - Mf [ - 3 (mJ - M 



m. 



Mf + MfMf-2 Mf 



m] 



(5Mf + 7Ml 



m„ 



J. 



abb 
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^ + ,ab 



-M} (Mf 



M} 



24 ml 



+ MfMf + 3(-ml + ml + Mf (Mf - Mf M, 



+ (Mf - Mff ( -2Mf + (ml + 7ml) Mf + (ml 
Mf (Mf 



2 I 2 
m^ + m^ 



m; 



l-ab 



M) 



4 ( m^ + m 



+ 5M? 



7Mf-2fm2 + 7 



m; 



m. 



7mg 



MJ - 2 (m^ 



5 M ? + 3 M' 



■QmlMf Mf (m| - 

+ (Mf - M|) (-3 m^ + M| + 2 Mf) 



m. 



m^ 



ml 



22 ml ml - 23 mf) Mf 
, 2" 



m, 



MfMf 



M-- 



l-ab 



abb ■ 
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mlMf 
p-q 



al3,^i,a/3 ^ -,ab 



-2Mj Mf 



Q((ml-Mf) 



mu 



Mf 
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+ 



+ (Mf - M|) (-8 Mt +{7 ml- 11 m^) + [ml - 



mi. 



2 ( 



5 m„ 



11 m^) Mj + 4 (m^-m^) j M; 



+ 2M| (^8M^- (m2-5m^)Mj-7 
- 6 Mj (^-3 + 2 (m^ + Mfj ml - 



m„ — mj 



12 Mj 
+ 4 M| 



Lab 



(m] - Ml 

J abb 1 



Mf 



{2ml- 2 



3 Mf) Mf 



48^^p(y ,Ft'^-' = M', 



{p ■ 



■2Mf-5{ml-2ml)M^ 



7mt 



8M| 



11 m^) ml + 4mi)Mi 



- ml (7 ml + ml) M} Mf - ml [ml - mff Mf 
+ [2 (mJ + {2ml - 5m^) 



Lab 



-2{2ml-mlml + mi)Mj + {ml 
-mlMf{7Mf -8{ml + ml) M 



+ (m 



ma m^ 
2 



Mf 



m\ 



Mf I 



+ 6mlMfMf {Mf -M 

J abb 1 



-2mlMf 



ml+ {ml-Mf + 2 Mf) m 



Lab 



48 



MP p(l) pa0,fi,af3 

- (llml + hmt\Mf 



-Mf 



m„ 



{Mf - Mf) (10 Mf 



m'l 



+ 



+ 3 {Mf -ml -3 ml) {-ml + ml + Mf) Mf 
2 (AMj -{5ml- ml) Mf + {ml - ml} 
+ Mf{ml + 10ml ml - llm^ - 5Mj 



Lab 



+ A{ml-2ml)Mf)Mf\h, 
+ 6 ml Mf Mf {ml + 3ml + Mf - 2 Mf) {Mf - Mf) J,,, , 



5ml - (Ami 



7M: 



m„ 



-ml + 2Mf + 5mlMf 



— 6 mf, 



Lab 



ml + ml + Mf) Mf 
{-2Mj+{ml-5ml)Mf 



m. 



m] 



Mf 



Lab 
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mlMf 



24 (ml 



24 {ml - mt) {ml -ml -2 Mf) {Mf - Mj) 

2 (m| -Mfy{AM^+{7ml + ml) + (m^ - ml 

2 (Mf - M| Mi)' (m/ - 2 (m^ + - 14 m^) M^ 



+ (m 



(m| - ) ' {mt Mj (mJ + 5 Ml 



M? 



3Mf) (m^ 



— m 



+ (m 



Mr 



4 ( 3 



3 — 2 m^ + 
3 M?! M? M? 



5 Mj) Mj + {9ml + 7Mj) Mf 



m„ 



- {Mf - M^y (3 {Mj -AMfMf + 3Mi 
+ 2 Mj (3 Mj - 4 Mf Mj + Mf + (4 Mf - 5 Mj 
+ M) (6 - 4 (mJ - 2 Mf ) + 7 Mj + 5 Mf - 12 Mj M/ 



+ Mj(l2 {Mf + Ml 
- 2 (A'f| - 2 Mf ) (mJ - Ml 



6 (5M) 



2Mf Mj 



(m 



m; 



3M^ 



(3M 



M- 



+ 4 M| (mJ + 8 Mf MJ - 4 Mf MJ + Mf 



mlM^f 



6{3Mf + 6MfMf 



2 ( M? + 26 M^ Mj - 31 Mf M| + 4 M^ 

,4 



Mf ) - (mJ - Mf ) {5Mf + Ml 



lab 



+ UmlM^ Mf [ - 2 (m^ - m^) Mf + (m^ - M^ 



+ (m 



m. 



m. 



m. 



{^l 
Ml 



2 m. — 3 m,. 



5 Mf ) (Mf - Mf) Mf 



M^ 



2 Mf + 2 + 2 mf - 4 m: M 



Mf - Mj)^ {ml ml + {2ml -3 K 



m^ + Mf 



J, 



abc ■ 



96 Mi Mi M P- 
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G 



12 (ml 



mi 



Mf 



M} 



2 ml Mf [12 {ml - ml) Ml 
-ml + ml + 2Mf)Mt 



+ (Mf - Mj M,)' (Mf - 2 (m^ + 9m2 - Uml) Mf + (m^ 
+ (mJ - Mf )' (4 Mf + (7m2 + m^) Mf + (m^ - rr 



m. 



M 



{Mf - 3 Mf ) {Mf - Mf ) ^ + {Mf - 3 Mf ) (mJ - Mf 
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(mJ - M^)^ ( (3 - 2 + 5 Mj) 



A [3 ml 



3 Mj ) M? 



/ 



9m2 + 7M|)M/ 



- (mJ - M^)^ (3 (m) - 4 M?^ M| + 3 
+ 2 Mj (3 - 4 Mj + M/ + (4 - 5 Mj 
+ (6 - 4 (mJ - 2Mf ) + 7 + 5 M/ - 12 
+ Mj (12 (m} + MI 



m„ 



6 (5M) 



2MfM} 



3M- 



-2 M| 



2M; 



3Mj 



+ ml Mj ( (mJ + 5 M^) (mJ - 



2 ( Af ? + 26 Mf Mi - 31 Af/ Af? 
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+ [m 
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A^/ 
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mi 



m^ 



(m] - Ml 
ml + 2ml-3ml-5M^) (M^ - M}) 



M 



Mt 



2Mf + 2[mi + 2mi-Amt]M, 



mi 



+ [m^ - [ml ml + {2ml -3 ^l) {^l + 



J, 



abc ■ 
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m. 



+ [Mf - Mj M,) (Mf -2(ml + 3ml - Sm^) + [ml 
MfY (-2 Mf + [ml + Imt) M^ + [ml - ml)' 



m. 



Mt 



2iM^f 



8 Aff Mj -AMfM^f+M^ 



m. 



+ [Mf [3 [3Mj + Q Mf M} - Mf) ml + [m] - M^f [m] + 4 Aff ) ) 



-ml [5M]-AMf) [M'f-Mf 



3 Af ? f 5 Af? 



7 



A^l 



M 



ml - M}f [m} - 2 Mf) [Mf - Ml 



+ 6MfM^ 



[Mf-M^) -2( 
+ ( mj, — mA ( + 2 — 3 m; 
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Qm^Mj [m^ + Mf) 



lab 



Mt 



M- 



Mf - M}) Mf 



Mf 



M- 



2 mt\ Mj 
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+ (ml - ml) (ml - ml) ) 



m. 



{^1 



2 ml 



3ml - Mf 



Mf 



J, 



ahc ■ 



48 M^M/P^y 

, 2 



mlM'j 



12 ( 



12 (m^-m^) Mf 
- Mj) Mf + (Mf - M] Mif (mi 
^nil)M^ + (ml-mfj^) 
-2M^+(ml + 7ml)M^ + (ml - ml 
+ mlMf[-2 (Mf + 8Mf Mj -AM^Mf + Ml 



3 ml 



(mI - Mj 



+ (m 



3{QMfM] 



Mf) ml + (mJ - M^)' (m| + 4 M; 



+ 9 



5 Mj - 4 



3mlMj (5 M 



(mJ - Ml 
M^y -Qm^^M* (Mf + Ml 



Mf - 2Mt ][Mi-Mi 



'aft 



+ 6 M; M/ 



+ [m 



mu 



m„ 



m. 



m. 



m. 



m. 



m„ 



Mf 



m. 



2 mi 



Mf 



m. 



(Mf - Mf 



M- 



3mi-M; 



2 { mi -2 m: 



Mf- 
Mf 



Mf 



Mf 



:) (ml 



mi 



2 ml -3ml- Ml 



(Mf - M' 



J, 



abc • 



(182) 



We turn to the contractions with the tensor P:^ of fll7ip . We obtain 



48' 



MiMf 



[p-q] 



3(2) 
af3,fi,al3 



A 



Mf 



-I2(m 



Mf) (ml + Mf) Mf 



12(Mf-Mf) (ml + ml + 3Mf)Mf 



.Q(-ml + ml + Mf) (Mf-M 



Mf 



m„ 



%(ml + Mf)ml + lml 



31 Mf + 22ml Mf 



l-ab 



m„ 



12 Mf 



24:mlMf 



mi + Mf) [mt + Mf)Mfh, 



m. 



Mf)Mf(Mf-Mf)jaab: 
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48 ■ 



(p-q) 



Mt 



12 (m^ - 



m. 



- 12 (M^ - Mj 



6 (m„ 



) [Ml 



Ml 



29 



10 m^M^ 



12 M| (M' 
24 m^M^ (m 

>(2) 



2 _ 



Lab 

4 



■nil 



M?) Mf 



7Mf + (7 ml -3ml -2 MJ) 



m 



t - lOm^ ml + 9m^ - 2 fm^ + m^) M?) 



4M; 



48 



/ 



12 mf 



Mf 



2 (m^ - llm^) + (ml - m^ 



+ 6 ml (m^ 



/ 



ml-ml + Mf 



lab 



+ 



{mi -2{ml + MJ) ml + {ml - Mf) ) Mf 
- Mj {ml -2{7ml + Mfj + 13 + M 
2AmlMj Mf {m] - Mf) J^,,, 



M? 



10 m^Mj 



48 



,(2) 



.B 



M]Mf 

ml - ml) M}) 
- 2 ( [ma - mbY 



a/3, /i, a/3 



Mf {-M^+{ml + ml + 3 M}) Mf 



m„ — 2 ( m; 



Ml 



rria + mbf - Ml 



Mf) ml + {ml - Mf 

4 + Mf) Mf 



lab 



lab : 



24 



M? 



,(2) 



5 



3M? 



+ (m^ + m^) Mj M^ - 2 (m^ 



3( 



mi 



Mf) Mf 



192 M^M^P^J .Cf'^bf 



12 -mlmA M] 



6 (m^ - m^) ml {ml - ml - Mf) {Mf - Mf) Mf 
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2ml{M]-Mfy{-2Mt+{ml-6 



1 1 
— m 



■2(m]- Mf 



Mf 



-Mt-(ml + ml)M^ + 2(ml- 
- (m^ - Mj)' ( - 3 Mf + (jml + + Mf 



2 

'b 



m^j [5 ml + 6 ml — m 
^^2Mj + 5M^M]-M^^..^^ 
+ [m} - Mf) (sMj-b Mf Mj + ml (s Mj - 2 M, 
+ Qm^^Mj - (ml + M'jml-2Mj) (^M] - Mf 



M?^ + {ml - m^ 

,4 



m. 



M] 



M] - Mf 



ibc 



+ 12 M| M/ [ {Mf 



ml — mlmh) 



ml -ml 



\2 



+ {ml - ml) (mJ - Mf 
{Mf - M\ 



' M?m 



J, 



abc 1 



192 MjMfP^l^^^^C^_%tf = Mf 
+ Q{ml- ml) ml {Mf - M? 



- 12 [ml - mlmt,) Mf 
{ml-ml + 7Mf)Mf 



ml - ml) M 



+ 2ml {M]-Mfy{-2Mf-{llml-Q 
+ {ml-ml)^)Mf 
+ 2 (Mj - Mf (-Mf - {ml + ml) Mf + 2 {ml 

- (Mf - M])' {Mf - {ml + Qml + 3ml) 
^ -ml){ml-6ml + 3ml)Mf+ ' 
(2 Mj + 5MfMf- Mf) ml 



m' 



— m 



Mf 

2 



+ 2m2Mf 



— t \ / ' 1 — / — t J 

- (mJ - Mf) (I6 M) + 5 Mf Mj + ml (2 Mf - 5 MJ 

+ 6 M^^ - {ml + Mj - 2 M^) (M^ - Mf ) ' 



mi. 



+ 12 M| Mf 



mlM'f 



3 M^ ((3 M^ + Mf ) + (m^ - 7 Mj) (mJ - M, 

2 - 2 Mf ) (mJ - Mf 
Mf)' 

Mj) (m^-m^ + 4Mf 



I- be 



-ml (-Ami 



■ mi. 



m: 



+ (mj, — 2m[ 



M) 



+ (ml-mimA M, 
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96 p. 



(2) 



c. 



a/3,^i,af3 ^ abc 



6 (ml 



m. 



m. 



3 m, 



-m^ + m; 
- {Mf-MfM,y{M^-{5ml 
-2(ml + m1+ [ml -3 ml 

2 



3ml + ml + M^) {m^ - Mj) 
2 + Mf) (m^ - Mf 



m. 



Mf 



5 M/ m; - 2 



,4 



(m; 

+ (mJ - Mf ) (-2 A'/^ + 5 M^Mj + ml{jM]-A M^ 



2mt 



Mfml 



Mf 



M) 



M 



mu 



+ 3mlMf (ml (m| + 3 M^) - (ml + Af|) (a/ J - Aff 
+ 6 Aff (m^ -ml + Mj - Aff ) ((m^ - m^) AfJ - (mft - A^f 
(^m^ + mi, m^) Mj — m^ {mi, + m,.) Mf 



J. 



48 



M\ 



f 



.(2) 



py^J £)Oil3,iJ,,al3 



abc 
, 2 



3mlM} 



-Mj + AMfMj 



Mt 



Af|-A^?^) (m^ + Aff)-2m^A^2(^AfJ + M; 



A^/ 



Af; 



12m^Aff + 2 m^ + Aff M 



(Af_ 

+ 12 ml Mj + (ml - Mfj ' (m} - Mf 



Mt 



lab 



+ Qml (ml - ml) M} Mf [Mf - (m^ + - 2 ml) M: 



-ml + mi + Mf] Af, 



48 



Af? 



,(2) 



.E. 



J abc 1 

2Mf 



(-Ami + 8^fe + Mfj M^ 



+ 2 (-m^ - 4 + 5 + (m^ + m^) Afj) Mf 



A^j(7m^ + 10 



a 



llml + 3Mf 



3 (3m2 



5 ml 



3 [ml -ml 



Mf (ml 



3 m, 



Mf 



-8Mf + (7ml - llml) Mf + (ml 



lab 



— mi 



Mf 



-Mf{- 7 ml + 5 [Mf -2mt] m^ + 17 ml + 2 Mf 



23 ml Mf ) M', 



Lab 



+ 6 ml Mf Mf (2 ml + 6 ml + Mf - 3Mf) (Mf - Mf) Jabb ■ 



m 



Mf) Ai", 



95 



96 m 



M 



f p(2) 



12 {{ml- 



mt]Mf 



3 {-ml + ml + Mf ) (Mf - M|) {-ml + ml + 7 M^) 
2 {m 



m; 



2[mi+ Ami) Mf + [m 



lab 



+ 



Mf 



(4 Mf - {5 ml- m^) + (m^ - n 
+ Mj - 2 (m^ + 19 ml) M| - 5 (m^ - 

+ M)(l6m^- (8m2 + 17M|) ml-Sm^ + M^ 



Mf 



37 mi Mf]M, 



lab 



-6 ml M] Mf (m| - Mf) [m/ - (-5 ml + 5ml + MJ) M?^ 



M; ( -9 + + 4 M; 



J, 



96 m 



P. 



(2) 



,ab 



M) 



24ml {-^l + ml + Mf) M, 



- (m^ - Mj) (m/ - 2 (m^ - 5m^) Mf + {ml - 
+ 2 (m| - ) Vm/ - 2 (m^ + 4 m^) Mf + {ml 



m, 



Mf 



m^ 



lab 



+ 



Mf 



{ - 2 Mf +{ml-5 ml) MJ + {ml - ml 
- Mj (-5 Mj + A{ml + 7 ml) M} + {ml - m^)') Mf 
+ 3Mf (-8 ml + 3M'jml- Mf + ml (s ml + Mj)) M/ 



lab 



+ 6mlM] Mf {m} - Ml 



M 



-3mi + 3mi + Mi) M- 



+ M] {3ml + 5ml - 2 M 



J, 



abb ' 



12 



,(2) 



E 



3 al3,fj.,al3 ab 



al3,fj,,al3 



— — 4 ml ml + 5 



M/ 



+ 2 m^-2m," M 



lab ■ 



p(2) pa^,^i,a/3 



3M? 



3 [m^ + m; 



m;)m 



48 



- (m^ + m2)MjMf + 2(m2 

m|M2 
[p-q] 



M^ 



p(2_) pal3,fj.,al3 
3 al3,^i,al3 +,ab 



3Mf-Mf)ml 



- (3 Aff + 7 Mj Mf + (mJ + 5 Mf )) + (9 Mf + 5 Mj M/ 
+ (l3Mf - 7Mj) + 2ml (s Mj Mf -23 M/ 



96 



- (ml - (5 + (7ml + uf) M^-5ml M 



12 ^(2) ^ 



(m^ - M]) (Mf + (m^ + ml 



2 (ml 



■ m'l 



QMfim 



ml 



96 m 



,(2) 



al3,ij,,al3 ab 



,af3,fi,al3 



- (m^ - Mfj (Mf - 58 ml M, 



' {p ■ q) 

- (?,ml-hSmlml-l7mf)Mf + 2(ml-?, 



M? 



(Ami- {2 



m. 



5m; 



+ 2mt))M^ 
2ml + Mf + 7mlMi 



+ (m 



(2ml + (l6m^ - 7Mf) + llm^M?^ 
14 + 28 + 5 M/) - 4 - 49 ml M, 



+ 12m2 (3m2 + ml - Mf) (-m^ + + M^) M, ^ 

[ - 2M) - 2 (2m^ - 9m^) + 8m^ (ml - 4m^) 
- 2 (m^ - 18 mlml + 9mlml + ^ 



ml) M} 

2 



Mf 



+ 



-(M; + 4(m^ + 4m^)Mj-5(m^ „^ ^ , 
+ (3Mj + 2 (4m2 + 5m2) - (l3m^ + 6mlml + 5m^) Mj 
+ 2(m;:-3m^m2 + 2m^))Mf 
QmlM]M^ (^Mf -M] 
3 ml + ml- M^f 



2 Ml 



-r4 



lab 

3 ml 



m 



4 m? 



J. 



a66 ) 



96 m 



+ 



2 p(2) pal3,ii,aP ^ 



24 



^2 MJ) 



+ (mJ - Mf)' (5M/ - [iml + 13m^) 

5 + 4 (m^ - 2 m() Mj + (m^ - m^)^ )Mf 



2 ( - 



+ Mj (13 Mf -2(7 ml + 13 m^) MJ + (ml 
2Mf (m^ - 14mg ml + ISm^ + iMf - 5 (m 



ml 



Mf 



ml\M))Ml 



lab 



+ 6mlMfMf (m1 -M 
'(5Mf + 3Ml 



3ml-M] -2 Mf) (m^ - M] 



■ mi. 



J, 



abb 



96 m 



2 p(2) pa^,IJ.,af3 ^ ^2 



a/3,ii,al3 -,o.b 



- 12 



m 



97 



+ Mf 



M] 



17 Mt + 2 {ml -5ml) Ml 



- [M'j-Mf 



m„ 



7M'f-A{2ml-ml)M} + (ml- 



13 ml) M^ 



7 {ml - ml. 



mi 



2 mi 



I 



mi 



Mf 



+ M; 5 - 2 (5 + 11 1 M; + 5 ( - m, 



Mf 



+ 6Mj {ml+ {2ml- M^j 
QmlM'jMf {m}-M\ 



m„ 



?,ml + ?>mlMj]M: 



2M. 



+ M) {5ml + 3ml - 3 M 



-ml + ml + M'j) M^ 



lab 
2 



abb • 



48 



{p ■ qY 



Mi 



AMf 



3 ml + 13 ml 



3Mf) M 



+ {9mt-mt + 3{mi + mt]Mj]M, 



— 2 m„ \ m„ — mt 



m2-ml + 3Mf 



lab 



+ Mf[2{- {2ml + ^ml) Mf + (m^ - 8 + 3 mt) Mj) 

+ 2 {ml- ma miy + (m) + 4 {ml + 4ml) M'j 



5 [m„ 



mi 



Mt 



lab 



+ 6 ml Mj Mf {Mf - Mj) {-ml + - Mj + 2 Mf) J^bb : 



48 



M]Mf 



[p ■ qY a0,fi,al3-^+,ab 



M] 



mf, + mf + Mf] M 



+ 2[Mj-Mf) (-M^-[mt + mt)Mf + 2[m, 
+ {Mf - Mf) {2ml + {^^l + ^Mf) ml - lOm^ - 7 Mf 
+ llmlMf)Mf 

- {-bMj + A{ml- 2ml) ^ 

- 2 Mj {Mj + {ml + ml) MJ - 2 {ml 



m. 



mi 



M? 



M 



lab 



+ 6mlMfMf{Mf-M]) {- 



mi. 



M) 



2Mf] J, 



'abb ■ 



A/f2 



Mj 



2Mf 



m„ 



17 ml + 



Mf) Mf 



llml + (10m^ + M;jm^ + m^ + 7m^M;jM^ 



-25m' 



4m2 ml 



m, 



u) 



lab 



+ 



6 fm^ +m2 - M^) MfM^ 



98 



[iMj-A (2ml - ^b) + (^a - ^ 



Mf 



- 6 ml M} [m] 



'-)(- 



-ml + ml- Mj + 2M^)j, 



ASmlM',P^l^^^G:!r' = 2mlM' 



f 



lab 

'abb 1 



— 6 [ml — m^ 



+ 



f 

{ml - ml) {Mf - Mf) (-ml + + 3 M^) Mf 



-2Mt+ [ml - 5m^) + (ml - ml 



'^mlmlM) {2M)^hMlM) - Ml ^ 
-2 ml ml (mJ - M^) (4 Mf + 5 + ml (2 Mf - 5 M 

(3ml + 2Mfml + 3Mj) (m] - M""^ 



'■ac 



- m 



c 2 4 
— m^ m^ 



+ 12 m^Mj + 



: m. 



3 Mj + Mf) Mf + 3ml (ml + Mf) (Mf - Mf 
' + (ml -Mff (ml + Mf) (Mf - Mf " 
(Mf - ml) (ml + 2 Mf) (Mf - Mf " 



+ (Mf -Mfyml + 6 m\ Mf (3 Mf 



-UmlM'Mfi- 



mu 



m,. 



+ (ml Mf - ml Mf) (- 



ASMlMlP^''^ ^Gff'^'"^ 



+ Mf - Mf 
-ml + ml + Mf 



(Mf - Mf 
ml -ml) (Mf 



lab 

Mf' 



m„ 



M: 



2mlMf 



2 2 
— m 



Jabc 1 

^Mf 



+ Mf 



+ 3 (ml - ml) (Mf - Mf) (-ml + ml + 3 Mf) Mf 
+ (Mf - Mff (-2Mt + (ml - Mf + (ml - 

ff [ (Mf - Mff ml + 3mt (ml + Mf) (Mf - Mf ' 
+ 2mlMf (2Mf + 5MfMf - M- 



m' 



4 c 4 2 
m,, — D m^ m^ 



(3M 



— 2 m 



(Mf - Mf) (4 Mf + 5MfMf + ml (2 M\ 
'^ + 2Mfml + 3Mf) (Mf-M, 



Mf) Mf 



bMf 



Mf 



+ 12 ml Mf + (ml - Mff (ml + Mf) (Mf - M 
+ 2 ml Mf (Mf - mf) (ml + 2Mf) (Mf - Mf) ' 
QmtMf (3Mf - ml) (Mf - Mf)' ' 
12 ml Mf Mf (ml - ml - 
+ (mlMf-mlMf)( 



■Mf + Mf 



lab 
/ 2 2 

[m^ - m^ 
Mf 



Mf 



Mf-M'f] mt 



abc 



24p15 .Hf^'^f" = M^f 



6 (ml 



m: 



' Mf 



3 (ml - mf) (-ml + ml + Mf) (Mf - Mf) M: 
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+ (m 



2{ml + M^)ml+{ml-M^ 



+ M^[(2Mf + 5M^ - M^) - (-5 + 9 + Mj) M; 



4\ 2 



6 



+ 6 -m^ M;M, 



Mi 



lab 

- Mf]m 



afec 



24 Mf P. 



(2) 



6 (m^ 



f al3,fj.,al3 abc "V^"/ 

Mi^f (rat- 2 H + + H " 



+ mlMf 



(2M^ + 5MfM| -Mf 
7 + 3 + M|) Mj + 2 (m^ + Mj) 



6 mt Mj - (ml - MfV (m} - ' 



l-ab 



+ 6 ml M} Ml [ (-ml + ml + ) (m^ - uf) (ml - 
+ M^ (ml - m,2)' + (ml - ml) ml (m| - 



J, 



abc • 



;i83) 
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